ApJ in press 



Heating and Ionization of X- Winds 

Hsien Shang^'^, Alfred E. Glassgold'^ , Frank H. Shu^ 
Astronomy Department, University of California, Berkeley, CA 94720 

o ■ 
o 

, hshangOcf a. harvard. edu,fshu, aeg@astro.berkeley.edu 

• and 

o 

'nI" ' Susana Lizano 

(N 



> 

cn 
in 
o 



Instituto de Astronomia, UNAM, 58090 Morelia, Michoacdn, Mexico 
lizanoSastrosmo . unam . mx 



ABSTRACT 

o '■ 

'q^, In order to compare the x-wind with observations, one needs to be able to cal- 

Q . culate its thermal and ionization properties. We formulate the physical basis for the 

streamline-by-streamline integration of the ionization and heat equations of the steady 
x-wind. In addition to the well-known processes associated with the interaction of stel- 
lar and accretion-funnel hot-spot radiation with the wind, we include X-ray heating and 
ionization, mechanical heating, and a revised calculation of ambipolar diffusion heating. 
The mechanical heating arises from fluctuations produced by star-disk interactions of 
the time dependent x-wind that are carried by the wind to large distances where they are 
dissipated in shocks, MHD waves, and turbulent cascades. We model the time-averaged 
heating by the scale-free volumetric heating rate, Tmech = apv^s~^, where p and v are 
the local mass density and wind speed, respectively, s is the distance from the origin, 
and a is a phenomenological constant. When we consider a partially-revealed but ac- 
tive young stellar object, we find that choosing a ~ 10^^ in our numerical calculations 
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produces temperatures and electron fractions that are high enough for the x-wind jet 
to radiate in the optical forbidden lines at the level and on the spatial scales that are 
observed. We also discuss a variety of applications of our thermal-chemical calculations 
that can lead to further observational checks of x-wind theory. 

Subject headings: X-winds, Jets, Herbig-Haro Objects, Young Stellar Objects, YSO 
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1. Introduction 

A refined and updated version of tlie disk-accretion paradigm for tlie formation of Sun-like 
stars has emerged in the last two decades through extensive observations and theoretical studies. 
In addition to the building up of the new star by accretion from a disk formed by the collapse of 
a rotating molecular cloud core, the formation of a low-mass star is accompanied by a remarkable 
bipolar OTitflow that can appear jet-like at optical and nearby wavelengths. Equally important is 
the crucial role of magnetic fields in retarding the initial collapse and in guiding both the accretion 
flow that feeds the star and the outflow that removes excess angular momentum. In addition to the 
strong evidence provided by the essentially universal detection of X-rays in low-mass young stellar 
objects (YSOs), magnetic fields have been measured directly with the Zeeman effect (e.g., Johns- 
KruU and Valenti, 2000). Although the general outline of a theory of low-mass star formation 
has emerged, the underlying mechanisms still need to be identified and understood, and strong 
efforts along these lines are in progress on a broad front, as can be seen in the reports at the recent 
conference Protostars and Planets /F (Mannings, Boss, and Russell 2000). 

One of the main goals of the theory is to develop a rational description of the active flows close 
to the central engine, i.e., the accretion funnel and the wind. Although there is a consensus that 
these flows are MHD in character, considerable disagreement exists over the specifics, as witnessed 
by the reviews of magnetocentrifugal winds at Protostars and Planets IV by Konigl and Pudritz 
(2000; disk winds) and by Shu, Najita, Shang, and Li (2000; x-wind), as well as the earlier review 
of wind theory by Pudritz and Rudcn (1993). Of course the only way to decide between alternative 
theories or to validate any particular theory is to make detailed comparisons with observations. 
Thus, it is the objective of the present paper to develop the basis for making such comparisons for 
the case of the x-wind (Shu et al. 1994a; Shu et al. 1994b; Najita & Shu 994; Ostriker & Shu 1995; 
Shu et al. 1995; henceforth Papers I-V). The x-wind model has the potential for understanding 
many aspects of low-mass star formation because it provides well-defined dynamical solutions that 
can be used to make detailed correlations and predictions of observational data. It should be 
clear that, because x-wind theory focuses on the inner region (or "central engine") of the star in 
formation of dimension 0.1 AU, observational tests require spectroscopy on the milli- arc-second 
spatial scale and better. Such observations are now becoming available with adaptive optics and 
interferometric techniques, as well as with the Hubble Space Telescope. 

As discussed by Shu et al. (2000), several crucial assumptions and implications of the x-wind 
model are supported by observations: the existence of a finite inner radius for a rapidly rotating 
inner disk; strong magnetization of the central star; magnetically channeled accretion; and phase 
relations between stellar rotation and the accretion funnel and the outflow. The model can also 
account for the large-scale kinematic properties of bipolar molecular outflows (Shu et al. 1991, 2000), 
and it has the potential to explain the optical observations of jets from young stellar objects. The 
latter possibility is a consequence of the remarkable property of the x-wind, whose 180° bipolar 
lobes self-coUimate towards the outflow axis into approximately cylindrical jets (Paper V). On this 
basis, Shang, Shu, &; Glassgold (1998, henceforth SSG) made synthetic images of optical jets that 
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bear a striking resemblance to the observed ones (see, e.g., Eisloffel et al. 2000). 

These conclusions have been made on the basis of the dynamical solution of the x-wind model 
which gives the density, velocity, and magnetic field configuration^. The images obtained by SSG 
were obtained by using constant values of the temperature and electron fraction based on previous 
analyses of the optical observations. The main goal of the present paper is to lay the foundation 
for obtaining the physical properties of the x-wind from first principles in order to calculate the 
distribution of the temperature, the ionization fraction, and other chemical abundances for the 
x-wind streamlines. These properties are required for calculating the fluxes of diagnostic lines and 
continua that observers can use to test the validity of the model. 

Even though the dynamics of the x-wind are largely decoupled from the thermal-chemical 
properties, the calculation of these properties for a 2-dimensional flow is very difficult. The closest 

previous work by Ruden, Glassgold, & Shu (1990, henceforth RGS; see also Glassgold, Mamon, &: 
Huggins 1991 for a parallel chemical study) done before the x-wind solutions were obtained, assumed 
spherical symmetry. They tried to anticipate some aspects of the 2-d axisymmetric solutions by 
modulating the radial density and velocity variations at small distances r. RGS found that radial 
winds quickly cool and become weakly ionized with increasing r. The importance of the run of 
ionization and temperature in outflows is well illustrated by the observations of optical jets, where 
phenomenological analyses of line strengths indicate that they are significantly ionized and hot, 
with electron fractions ~ 0.01 — 0.1 and temperatures T 5,000 — 10, 000 K (e.g., Bacciotti 
2001). Shocks have long been the favored mechanism for producing these conditions (e.g., Raga, 
Bohm, and Canto 1996; Hartigan, Bally, Reipurth, and Morse 2000). Aside from the "final" bow 
shocks with the ambient medium, associated with the Herbig-Haro objects, it has been unclear 
how the central young stellar object (YSO) can affect the global properties of the jet at large 
distances from the source. Bacciotti et al. (1995) suggested that the jet retains a high level of 
ionization characteristic of the source region by virtue of slow radiative recombination in the flow. 
Our calculations support this idea, once we add an important missing ingredient, the ionization 
of the base of the wind by the X-rays that are observed to accompany essentially all YSOs (e.g., 
Feigelson Sz Montmerle 1999). 

Not only are the X-rays effective in ionizing the wind close to the source, they also help maintain 
the ionization level at large distances. The emissivity of the flow is determined by the temperature, 
and it is essential for understanding jets to also be able to achieve the high temperatures indicated 

by the observations of the optical forbidden lines. Here shocks can play an important role, and we 
will develop a global model of wind heating based on stochastic shock dissipation. Although our 
model of mechanical heating is supported by MHD simulations (e.g., Ostriker et al. 1999), at this 
stage it is essentially phenomenological. By adopting physically reasonable parameters, we will be 
able to obtain images of jets that resemble the observations. We should also be able to relieve the 
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difficulties encountered by RGS, since the wide-angle component of the x-wind seems Ukely to be 
warm and moderately ionized. In a reconsideration of ambipolar diffusion heating, we will find that 
a new atomic coefficient provides a much reduced role for this process in warm atomic regions, in 
disagreement with earlier results by Safier (1993). 

The main goal of this paper is to develop a coherent thermal-chemical foundation for the x- 
wind. We build on the previous study by RGS, but add important new processes for ionization 
(X-rays) and heating (mechanical or turbulent shock heating) . As the main illustrative application, 
we consider the jets as observed in optical forbidden lines. The rest of this paper is organized as 
follows. In the next section §2, we give a general description of the model, and in the following 
sections we focus on X-ray ionization (§3), ambipolar diffusion heating (§4), and mechanical heating 
(§5). Modeling results are then presented in §6 for the case of an active solar-mass YSO in a partially 
revealed phase. We then discuss further implications of the calculations in §7 in the context of future 
detailed studies that bear directly on observations. The paper contains a number of appendices 
that supplement the technical basis for §§2-5. 



2. Formulation of the Model 

The calculation of the thermochemical properties of the x-wind is separable from the dynamical 
problem (treated in Papers I-V) in the cold ideal MHD limit. As long as the electron fraction is 
large enough and the thermal energy is small enough (for the thermal pressure to be small compared 
with the kinetic and magnetic energies), the MHD approximation can be made. We of course check 
that the thermochemical calculations reported here satisfy these assumptions. Most of this section 
will be devoted to formulating the thermal and chemical equations that need to be solved along 
each streamline. As discussed in §2.1, the streamlines are obtained on the basis of the dynamical 
solution for the x-wind obtained by Shang (1998). 



2.1. Dynamics 

The exact numerical solution for the x-wind in Paper III does not provide a practical basis 
for thermochemical modeling because it is restricted to the sub-Alfvenic region. Similarly, the 
asymptotic solution in Paper V does not apply at small distances. We use a semi-analytic approach 
developed by Shang (1998), where a global x-wind solution is obtained by interpolating between the 
two extreme solutions developed in Paper II and Paper V. This solution applies to the steady and 
axisymmetric X-wind flow that we model here. It satisfies the conservation laws of mass, specific 
angular momentum, and energy, which are expressed in terms of the conserved quantities 
J{ip), and H{il)) on the streamlines if) as shown in Paper II. In steady state, the field lines co-rotate 
with the star at angular velocity VL.^,. Unlike Papers II and III, which employ a reference frame that 
rotates with the stellar angular velocity il*, we describe the flow in an inertial frame. We follow 
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Paper I in using non-dimensionalized equations based on the units for length, velocity, density, and 
magnetic field: R^, i^xRx, Myj/ATrR^^l^, {VL^M^ / R^)^/'^ , where Q.^ = is the angular velocity 
at the X-point R^, and is the mass- loss rate of the x-wind. Our calculation is restricted to 
the case treated in Papers I-V of aligned magnetic and rotation axes. Many of the same physical 
processes discussed in the following sections would be operative in the more general case, but the 
required dynamical solutions do not exist yet. 

In the asymptotic regime, an x-wind streamline ip is represented in spherical coordinates r — 9 
by the parametric equations, 

r„ = ^cosh[F(C,l)] (2-1) 

and 

0„ = sin-i{sech[F(C,V')]}, (2-2) 

where 

and the subscript a stands for asymptotic. The stream function ijj has the range < ■0 < 1 and 
labels the mass fraction carried by the x-wind from the horizontal plane. The quantity /3 is the 
ratio of the magnetic-to-mass-flux, i.e., MHD field-freezing in the co-rotating frame implies that the 
magnetic field and mass flux are proportional to one another: B = f3pu. The condition V • B = 
requires that /? be conserved on streamlines, i.e., P = Similarly, the conservation of total 

specific angular momentum requires that the amount carried by matter in the inertial frame, wv^p, 
plus the amount carried by Maxwell torques, — zi7i?^B, per unit mass flux, p\i, is a function of V' 
alone: 'UJ[v^p — (3"^ p{v^p — w)\ = J{ip). The function C is a slowly decreasing function of r introduced in 
Paper V that vanishes logarithmically when r ^ oo, consistent with the vanishing of the current at 
infinity, independent of streamline. The locus ijj = 1 is obtained by applying approximate pressure 
balance, 3^ = 3"^, across the wind-deadzone interface. 

In the cold limit, the function f3{'ijj) cannot be chosen completely arbitrarily; otherwise the 
magnetic field, mass flux, and mass density will diverge on the uppermost streamline {ip 1) 
as the x-wind leaves the x-region. For modeling purposes, Shang (1998) adopts the following 
distribution of magnetic field to mass flux: 

/3(V) = /3o(l-^)"'/', (2-4) 
where /3o is a numerical constant related to the mean value of (3 averaged over streamlines: 

P = jy{i;)d^ = ^Po. (2-5) 

The singularity in equation (2-4) is of no real concern: it reflects the fact that the magnetic field 
B = /3pu is nonzero on the uppermost x-wind streamline, where by definition p must become 
vanishingly small while u remains finite. 
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In order to make contact with the inner solution, we follow Paper II and use pseudopolar 
coordinates (s, (p, i")) with the origin of s at the x-point and the angle i? measured in the meridional 
plane starting from zero at the equator: 

m = 1 + s cos z = s sin -d. (2-6) 

When r 1, s r and ■& — 7r/2 — 9. Equations (2-1) and (2-2) now become, 

s„ = ^cosh[F(C,l)], (2-7) 
i?a = cos-^{sech[F(C, ^)]}. (2-8) 

Near the x-point, the magnetic field lines emerge uniformly and define a fan within a 60° sector 
above the equatorial plane (0 < t?/ < 7r/3), 

^/(V)=^/(0)i r/3(^W, (2-9) 
P Jo 

where the subscript / stands for fan. The uppermost streamline with tp = 1 emerges at "df = 7r/3 
and the lower-most streamline with ^ = emerges at -iJj = 0. 

The interpolation between the outer asymptotic and inner fan solutions is based on the equa- 
tions, 

siC) = ^ cosh[F(C, 1)] - ^ cosh[F(C„ 1)], (2-10) 

and 

^C, ^) = a^fiiP) + (1 - a)MC, (2-11) 

The second term in equation (2-10) (with the constant Cx) has been introduced so that the function 
s = s{C) reduces to zero at the x-point {s{Cx) = 0). The reference value Cx is chosen so that 
equation (2-2) yields 7r/3, i.e. , the asymptotic streamlines occupy the same angular region as the 
fan emerging from the x-point (and are in approximate pressure balance). An interpolation function 
a(C, i)) was chosen to make a smooth transition between the fan and the asymptotic solutions. 



2.2. Heating and Ionization 

The equations governing the temperature T and electron fraction are 

-— = T— lnn(l + Xe + XHe)"^/^ + = (6? - >C) , (2-12) 

2 Dt Dt ^ e Hey 1 + + ^jj^ y > y J 

^ = V-V. (2-13) 
Dt ^ ' 

where = v • V is the substantial time derivative along a streamline, V = P/n and V = D/n 

are ionization production and destruction rates (dimensions s~^), and Q = T/nk and JC = A/nk 
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axe heating and cooling rates (dimensions Ks~^). The Latin symbols P, D, G, and L are the 
corresponding volumetric rates, and n is the total number of hydrogen nuclei per unit volume 
(usually denoted nn in the literature). We use n to define abundances, e.g., = ric/n is the 
electron fraction, XHe is the He abundance, and (1 + Xe + X}ie)i^ is the total number density of 
particles ignoring elements heavier than He. Many of the source terms in equations (2-12) and 
(2-13) were formulated in a unified way by RGS, and we use their expressions for individual terms 
in P, D, G, and L whenever possible. We discuss those terms in RGS which require change in §2.4 
and appendix A and B, and new ionization and heating sources in §§2-5. Table 1 lists the processes 
included in the present study. 

Two basic time scales are those for adiabatic cooling and recombination. Near the source, 
their values for the fiducial model introduced in §6 are of the order of tens of seconds and days, 
respectively. We usually omit explicit consideration of the effects of the ionization of He, and adopt 
a^He = 0.1 in calculations described in this paper. In the regions of main concern, the central part of 
the wind in and around the jet, the dominant ionizing agents are the secondary electrons produced 
following X-ray ionization of heavy elements. As long as Xe > 3 x 10~^, the inclusion of He would 
mcrease Xe by the factor (1 -|- XHe)^/^ or 5%. 



2.3. The Accretion Hot Spot 

The x-wind may be heated and ionized by several sources of radiation. In addition to the 

photosphere, these include a "hot spot" produced by the infall of the accretion funnel onto the stellar 
surface; the accretion disk, which scatters and re-emits absorbed stellar radiation and radiates 
energy generated by viscous dissipation; and the star-disk magnetosphere, which emits thermal 
and nonthermal radiation in both soft (< IkeV) and hard (> IkeV) X-rays. Because we are 
mainly concerned with the upper streamlines (close to the axis) that constitute the bulk of the 
mass of the inner wind, we can safely ignore the disk radiation. Not only is the mean photon 
energy smaller than for stellar radiation, the radiation emitted by the disk is significantly diluted 
by geometry and absorbed by the wind. The effects of disk radiation and gas-dynamic interactions 
with the disk atmosphere are important for the lower streamlines, which constitute the wide angle 
wind. We focus here on the hot spot radiation and treat the X-rays in detail in §3. 

We model the hot spot following Ostriker &i Shu (1995). For an axisymmetric x-wind, the hot 
spot is an annulus located between co-latitudes 9i and 62. The hot spot luminosity is 

i?? sin^ 3R, 



-^h = (1 - /) 5 

-ft* 



1 + 



2El 2R, 



(2-14) 



where / is the fraction of the disk accretion rate Md that is lost by the x-wind (and thus (1 — /) 
is the mass transfer rate of the funnel flow). The hot spot covers a fraction = cos^i — cos02 of 
the surface area of the star, and the mean colatitude of the hot spot is sin^h = Fh/{02 — Oi)- We 
use parameters for the fiducial case in §6 that are appropriate for the preferred model of Ostriker 



Table 1. Physical Processes 



Production of Ionization 



H~ photodetachment §2.4 

Balmer continuum photoionization §2.4, App. B 

Electronic coUisional ionization §2.4 

X-Rays §3, App. C 



Destruction of Ionization 
Radiative recombination §2.4 



Heating 

Photodetachment of H~ §2.4 
Balmer continuum photoionization of H §2.4 

H"'"-H~ neutralization §2.4, App. A 

Ambipolar Diffusion §4, App. D, App. E 

X-Rays §3 

Mechanical §5 



Cooling 

Adiabatic §2.2 

H~ radiative attachment §2.4 

Recombination of H"*" §2.4 

Lyman a §2.4 

CoUisional ionization §2.4 

Heavy element line radiation §2.5 
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& Shu: 9i = 26.6° and 62 = 33.2°, and with R^/R^ = 0.20, Fy, ^ 0.06 and sin^ 6*1, ^ 0.5. Calvet & 
Gullbring (1998) and Gullbring et al. (2000) have analyzed the spectral energy distributions of T 
Tauri stars with an accretion shock model based on a dipolar magnetic field. They find that lightly 
veiled T Tauri stars have Fh in the range 0.001 — 0.05. Such small hot spots can be realized with 
a generalization of x-wind theory that brings in higher multipole fields (Mohanty & Shu 2001). 

The presence of the hot spot is a potential complication in the radiation transfer needed 
to calculate the absorption of radiation by the wind. We could achieve some simplification by 
exploiting the axial symmetry of the steady x-wind and replacing the hot spot by equivalent sources 
located at the origin and along the polar axis of symmetry. Because we are mainly interested in 
radial distances r » i?*, however, wc can simply replace the hot spot by a source located at the 
origin and calculate the mean intensity of the hot spot radiation as a blackbody of temperature Th 
with the approximate dilution factor. 



W^ = F^^^. (2-15) 



The effective temperature Th is given by 



T^ = Tt+[ ^Top (2-16) 



KF^A-KRla 

if we assume that the accretion shock structure is optically thin to the underlying stellar (photo- 
spheric) continuum and denote the Stefan-Boltzmann constant by ctb- The corresponding dilution 
factor for the stellar photospheric radiation is 

^* = (l-^0^- (2-17) 

The mean intensities of the stellar and the hot spot radiation fields before they enter the wind are 
obtained by multiplying the dilutions factors by ctbT^ and (J^BTf, respectively. Although neither the 
stellar nor the hot-spot radiation is accurately represented by a blackbody spectrum at wavelengths 
shortward of the Balmer continuum, the blackbody approximation does not lead to any serious error 
in this paper because the X-rays are the dominant ionization source (§3). 



2.4. Radiative and CoUisional Processes for Hydrogen 

The stellar and hot spot radiation fields are important in the inner wind for detaching the 
electron of the negative hydrogen ion (in the continuum shortward of A = 1.647 //m) and for pho- 
toionizing the n = 2 level of atomic hydrogen (in the Balmer continuum shortward of A = 0.365 jJLva). 
These processes and their inverses were discussed in detail by RGS for a stellar blackbody, and we 
describe only the most important changes in their methodology. RGS expressed the rates for these 
photoionization processes in the form, 

r = WgAX, (2-18) 
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where x is the abundance of the species (H^ or H(2), where we write H(n) for the H atom with 
principal quantum number n), is the standard dilution factor for a star, and gp^ is the absorption 
rate for the appropriate continuum. Applying the results of §2.3, our ionization rates are, 

Vk = {W, gA,* + VFh 9A,h) X, (2-19) 

where VF* and Wh were defined in §2.3 and gA,* and ^A.h arc absorption rates evaluated using 
blackbody spectra at the stellar and hot spot temperatures, respectively. 



2.4- 1- The Negative Hydrogen Ion 
The rate equation for H~ is: 

-^x(H-) = xe(l - xMH-)n - [Wi,gA,h{H~) + W^gA^H') + x^nk^] x(H-). (2-20) 

In addition to photodctachnicnt, is also destroyed rapidly by neutralization with H+ in which a 
hydrogen atom is produced in an excited state, usually n = 3 for energies less than 1 eV (e.g., Fussen 
k, Kubach 1986). The inverse reaction, H(l) + H(3) H^+H+, contributes little to the production 
of H~, mainly because the abundance of excited H(3) atoms is so small. The potentially important 
reaction, H(l) -|- H(l) H~ -|- H+, is endothermic by 12.85 eV and has a relatively small cross 
section even well above threshold. It can therefore be safely ignored as a production mechanism 
for H^. We use a rate coefficient k± somewhat larger than in RGS, as explained in Appendix 
A. We calculate the photo rates in equation (2-20) by numerically integrating the Wishart (1979) 
photodetachment cross section a,y(H~) in RGS Eq. (A16), 

and in RGS Eq. (A18), 

gA,.(H",r.) = At, ""^j^ ^ B4T,)du, 5A,h(H-, Th) = 47r ^ "'^^^ ^ B4Th)du, (2-22) 

where is the threshold frequency for H~ photodetachment at A_ = 1.647 /xm. The subscript 
S in equation (2-21) labels as a spontaneous rate. We use detailed balance (RGS Eq. [20]) to 
obtain the rate coefficient for photodetachment from , 

fcs(H-,r)^^e'^-°/'=^3s(r). (2-23) 

We note that g'A,*(H~, T*) and yA,h(H~, Th) are functions of the stellar and hot spot temperatures, 
whereas gs(il~,T) is a function of kinetic temperature. We need to calculate 5A,*(H~,r*) and 
5A,h(H'~, Th) just once at the beginning of a calculation, but we need to know A;s(H^, T) everywhere 
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on each streamline. In order to save computing time, we have fit ks by a formula similar to that 
introduced by Stancil, Lepp, & Dalgarno (1998), 

A;s(H-,T) = 1.33 x IQ-^^ (T^-^^ + r0-4)(e-r/9320 + ^-T/isooO)^ ^2-24) 

which reproduces the numerical calculations to better than 15% over the temperature range from 
10^ — 2 X 10^ K. We have ignored attenuation in equation (2-22) because the optical depth in the 
H~ continuum is small due to the low abundance of H~ . 

Associated with the gain and loss terms in the rate equation for H~, equation (2-20), are the 
heating and cooling rates, 

g{R-) = [T^h5A,h(H-,Th)TA,h(H-) + W*5a,*(H-, r,)TA,*(H-)] x(H-), (2-25) 

£(H-) =Xe(l-Xe)nA;s(H-,r)rs(H-,r). (2-26) 

The heating rates (in K units), (/A,h(H~, Th)TA,h(H~) and yA,*(H^, r*)rA,*(H^), are given by inte- 
grals like those in equation (2-22) except that an additional factor {y—uJ) appears in the integrands, 
as is appropriate to photoelectron heating. The cooling rate fcs (H~ , T)T^ (H~ , T) is given by a sim- 
ilar modification of the equation for go,. Again wc evaluate TA,h and Ta,* numerically once per 
calculational case and use the following fit for Ts (H~ , T) 

Ts{ii-,T) = T (l^^l^) (e-^/™ + e-^/^°°°) , (2-27) 

where ^- = T^/T = 8750 K/T; this approximation is accurate to within a few percent over the 
temperature range from 10^ — 10^ K. 



2.4-2. The n = 2 Level of Atomic Hydrogen 

Our treatment of the H atom is essentially standard Case B recombination theory modified for 
a 2-d axisymmetric wind by use of the Sobolev approximation (e.g., Shu 1991, Chapter 9). It goes 
beyond RGS by including electronic collisional and X-ray ionization and excitation. The broader 
implications of the X-rays are taken up in §3. We adopt a simplified two-level plus continuum 
model for the H atom, and formulate and solve the steady-state population equations in Appendix 
B. The population of the n = 2 level is given in the convention defined after equation (2-12) by, 

(1 - Xe) „ k2inXc + ^21/321 

l + Q ki2nxc + kicUXc + Cic + C12 

where A21 is the spontaneous radiative decay rate of the Lyman-a line, /32i is the corresponding 
escape probability, 

1 — 

P21 = , (2-29) 

T21 
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and T21 is the (locally-calculated) optical depth of the line, 

T2i = A2i/i2(l - a^ejn-S' , S=- ^o^- (2-30) 

The rate coefHcients ki2 and k2i describe electronic coUisional excitation and de-excitation of the 
n = 2 level; kic is the rate coefficient for electronic collisional ionization of the n = 1 level. The 
quantities C.ic and (^12 axe respectively the rates at which X-rays ionize and excite the H-atom ground 
state, as described in §3 and Appendix B. The main additional changes from RGS are that the 
radial coordinate r is replaced by the streamline coordinate s and the characteristic length 2f„,/3r 
is replaced by S in equation (2-30). 

Our treatment of photoionization by the Balmer continua and the recombination into levels 
n > 2 closely parallels the above formulation for H~. It is actually simpler because the integrals over 
frequency can be done in closed form when the cross section has the Kramers dependence (RGS 
Appendix B). The coefficients in the photoionization rates for the stellar and hot spot blackbodies 
in equation (2-19) are then given by integrals like those in equation (2-22), with ai,{'H.~) replaced 
by the Kramers cross section (RGS Eq. [21]); the result is 

00 00 
5A,/»(2) = e-^''^^2ulai ^ Ei{jT2/T^), qaA'^) = e-^'^^2ulai ^ E,{jT2/T,), (2-31) 

i=i 3=1 

where U2 is the frcc|iiency threshold for photoionization from the n = 2 level and T2 is the same 
quantity in temperature units. The factor exp(— T2c) reflects the fact that we include absorption of 
the Balmer continuum, which is well calculated by replacing the frequency dependent attenuation 
factor exp(— r,^) by its threshold value. The optical depth at threshold is 

T2c = 2aiN2, (2-32) 

where ^ 

N2= X2ndr, (2-33) 
Jo 

is the radial column density for the n = 2 level. 

The rate at which radiative recombination occurs is given by, 

Ve = xlnks{>2). (2-34) 

We calculate the rate coefficient ks{> 2) from RGS (their Eq. [B15] and preceding equations). As 
discussed by RGS, the evaluation involves replacing sums over the principal quantum number n from 
n = 2 to nmax by integrals. RGS determined nmax from the condition derived by Seaton (1964) 
for the equality of the rates for collisional ionization and spontaneous decay for large quantum 
numbers. The Seaton condition characterizes an important property of the large n level population 
(where the rate of increase of the departure coefficients begins to decrease), but it does not describe 
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the size of the atom, i.e., rimax- In fact, for the conditions in the wind, nmax is determined by Stark 
broadening (e.g. Mihalas 1970), 

nms« = 1234 (ne cm^)-^/!^. (2-35) 

Despite the small exponent, nmax ranges from ~ 30 at the base of the wind to 100 — 200 at large 
distances. This variation has little effect on the value of the recombination rate coefficient, but it 
does have a bearing on the observability of sub-mm radio recombination lines (§6). 

The photoionization of the n = 2 level and the recombination to levels n > 2 lead to heating 
and cooling, respectively: 

g{2) = [W^h5A,h(2)TA,h(2) + W*5A,*(2)rA,*(2)] X2, (2-36) 

/:(> 2) = x^nA:s(> 2)Ts(> 2). (2-37) 

The heating energies, Ta,*(2) and rA,h(2), are calculated from the closed form expression, RGS 
Eq. (B8), and the recombination cooling A:s(> 2)Ts(> 2) from RGS Eqs. (27), (B13), and (B16). 

The heating and cooling by the H atom occurs by collisional excitation, de-excitation, and 
ionization and by photoionization and electronic recombination. When the net effect of collisional 
excitation and de-excitation are expressed in terms of the standard formula for Lyman-a line 
cooling, 

jO-Ly-a = X2A21P21T21. (2-38) 

the heating is given by (Eq. [B12]), 

g-C = -CLy-a + G{2) - C{> 2) - C^ow + Gk, (2-39) 

where 

-CcoU = ^XeikicUi + k2cn2)Tic = ^Xe{l - Xe)nTic + Y^I^) ' C^-^O) 

and 

= 2.22Cri2; (2-41) 

Tic and Tu are the ionization potential and the excitation energy of the first excited level of atomic 
hydrogen (in K units) , and kic and k2c are the rate coefficients for collisional ionization of the n = 1 
and n = 2 levels. As discussed in Appendix B, equation (2-41) defines an indirect X-ray heating 
term that arises because we use equation (2-38) to eliminate the net heating from collisional heating 
and cooling in favor of Lyman-a cooling. This term arises from the effects of the X-rays on the 
population of the n = 2 level of atomic hydrogen. The collisional ionization rates in equation (2-40) 
are discussed in Appendix B. 



2.5. Heavy Element Cooling 



For the temperatures important in the inner wind, the heavy elements mainly contribute to 
the cooling by forbidden line transitions, many of which serve as observational diagnostics for the 
jets that emanate from YSOs. We have included the forbidden transitions of 01, SII, and NI at 
solar abundances using 5-level model ions as in SSG. N I cooling is about one order of magnitude 
smaller than 1 and S 11 cooling in our calculations. We have not included Fe 11 which, by virtue 
of a large number of transitions, can be an important coolant at high densities. Only in the last 
few years have realistic calculations of Fe cooling been published (Woitke and Sedlmayr 1999; 
Verner et al. 1999, 2000). According to the cooling functions displayed graphically by Woitke and 
Sedlmayr, Fell cooling should be small compared to the dominant adiabatic cooling in our wind 
models. The reason for this is that, for the densities characteristic of the x-wind jet, the strength 
of the Fe 11 cooling is determined by the abundance of iron compared to oxygen. It is only at much 
higher densities, where the Fc 11 lines become optically thick, that the very large number of Fe 11 
lines makes this ion a more powerful coolant than more abundant ions with simpler level structure. 



3. X-Rays 

The ionization and heating of the x-wind by X-rays are another example of how X-rays pro- 
duced by YSOs affect the physical conditions for star formation. Earlier we used observations of 
YSO X-rays to determine the flux of energetic particles that can produce short-lived radionuclides 
at the time the Sun was formed (Shu et al. 1997; Lee et al. 1998; Gounelle et al. 2001). We have also 
found that the observed level of X-rays is more than sufficient to provide the dominant source of 
ionization in the atmospheres of protoplanetary disks (Glassgold, Najita, &: Igea 1997, henceforth 
GNl; Igca & Glassgold 1999). Here we show how these X-rays affect the degree of ionization at the 
base of the jets observed in low-mass stars. 

The ROSAT measurements of nearby clusters of newly formed stars provide us with a good 
idea of the average soft X-ray flux from T Tauri stars, and ASGA has extended these results to 
higher X-ray energies albeit at lower angular resolution (e.g., Feigelson & Montmerle 1999). The 
mean soft X-ray luminosity for a typical T-Tauri star is Lx ~ 5 x 10^^ erg s~^. The satellite 
observatories have also begun to elucidate the dependence of the X-ray emission on the age of the 
YSO, the amount of variability, and the properties of their flares. Although soft X-ray observations 
are often fit by a single thermal spectrum, there arc good reasons to expect that the X-ray spectra of 
YSOs are more complex. This is certainly the case for (magnetically active) M-dwarf stars (Schmitt 
et al. 1990). Some high-quality ROSAT spectra of YSOs also require a two-temperature fit (e.g.. 
Skinner & Walter 1998). The ASGA YSO spectra also indicate that there is generally a hard as 
well as a soft component. The simulations of the Sun as an X-ray star (Peres et al. 2000), based on 
modeling of YOHKOH Soft X-Ray Telescope observations, suggest that the high-energy component 
seen in low-resolution astronomical satellite measurements of YSOs is associated with small and 
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large flares. This situation is quite understandable since hard X-rays are released by energetic 

flares generated by reconnection events, whereas soft X-rays are mainly emitted by coronal gas 
trapped by large loops of closed stellar magnetic field lines. Another conclusion indicated by 
the available data is that the younger YSOs are stronger and harder X-ray emitters than the 
older ones. Again this is in accord with theoretical ideas about the effects of accretion on the 
magnetospheres of YSOs, which are prone to increased magnetic activity during epochs of enhanced 
accretion. Of course many of the limitations in our information on YSO X-rays will be removed 
by measurements of increased sensitivity and angular resolution now being made with the Chandra 
X-ray Observatory and XMM. The Chandra observations (Garmire et al. 2000; Fcigelson et al. 2001, 
private communication), already support the conclusion from ASCA (Koyama 1999; Tsuboi 1999) 
that there are significant differences between younger and older YSOs with respect to the hardness 
of spectra and the frequency of flares. 

In modeling the X-ray spectra of YSOs, we assume that a soft and a hard X-ray component are 
present with temperatures kTx ~ 1 keV and kTx ~ 2 — 5 keV, respectively. Since we focus on the 
more active YSOs that produce optical jets (and correspond to infrared classes I and II), we assume 
that the luminosities of the two components are the same order of magnitude, and that the total 
X-ray luminosity is substantially larger than the typical revealed Tauri star, e.g., Lx ~ 10^^ erg s~^. 
Of course Lx can be even higher during large flares (Feigelson & Montmerle 1999; Tsuboi 1999; 
Stelzer et al. 2000). As in the discussion of the hot spot radiation in §2.3, the calculation of the 
effects of the X-rays is complicated by the extended nature of the emitting regions and by the lack 
of relevant observational information on their physical properties. According to the x-wind model 
(see Fig. 1 of Shu et al. 1997), soft X-rays can be expected to arise within the region of closed 
stellar magnetic field lines (of linear dimension several i?*) and from the region above the funnel 
flow and beneath the helmet dome and kink point (of linear dimension .Rx)- Hard as well as soft 
X-rays may be produced by reconnection events that occur either in the reconnection ring (in the 
equatorial plane from 0.75 to i?x) or along the helmet streamer above the kink point (at radial 
distances ~ i?x)- We deal with this complicated geometrical situation by representing the several 
finite sources by a set of axial point sources: half of both the soft and the hard X-rays are assumed 
to emanate from the origin and the other half from points displaced along the z-axis by ibl.Oi?x- 

We follow the calculation by GNI of the X-ray ionization rate, which is based on the energy- 
smoothed cosmic photoelectric absorption cross section per H nucleus based on the compilation of 
Henke, Gullikson, and Davis (1993), similar to that used by Morrison and McCammon (1983), 

ape{E) = (7pe(l keV){keV/Ef; (3-1) 

for solar abundances, p = 2.485 and (Tpe(keV) = 2.27 x 10~^^ cm^. We assume that moderate 
and high-energy X-rays are most important (E > 1 keV), and we ignore the relatively small 
contributions of the primary and Auger electrons compared to the dominant secondary electrons. 
We also assume that the electron fraction and the temperature are not high enough for the heavy 
atoms to be very ionized, so that the main X-rays absorbers are the K- and L-shells of heavy 
atoms. In other words, we approximate the production rate for primary photoelectrons as if the 
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heavy atoms are in their ground states using equation (3-1). In the format of equation (2-13), the 
ionization rate due to the secondary electrons is 

where Lx{E) is the X-ray luminosity per unit energy interval, eion is the energy to make an ion 
pair (about 36 eV for an unionized hydrogen-helium mixture, according to Dalgarno, Yan, &; Liu 
1999), and rx is the X-ray optical depth. 

rr 

Tx = apeikTx)N, N= ndr, (3-3) 

Jo 

We introduce a low-energy cutoff Eq because the smoothing used to obtain the power-law fit, 
equation (3-1), removes the thresholds in the underlying photoelectric cross sections. The smallest 
threshold is 0.0136 keV, but the operative cutoff may well be larger due to absorption in the source. 
Defining ^ = E/kTx, we choose = 0.01, e.g., Eq = 0.1 keV for a thermal X-ray spectrum with 
kTx = 1 keV because. Below this energy, the secondary electrons no longer dominate the ionization. 
We also need the X-ray heating rate, 

Gx = ^J^ ^^^<7pe{E) e-^ {y^,,tE) dE , (3-4) 

where yheat is the fraction of the X-ray energy that heats the gas. Both eion and yheat are functions of 
energy E and electron fraction Xe, the latter because of Coulomb scattering between the secondary 
and ambient electrons. We exploit the fact that, at energies much larger than characteristic atomic 
energies (measured by the ionization potential), ejon and t/hcat arc approximately independent of 
energy. We can then extract the asymptotic factor l/eion from equation (3-2) and the asymptotic 
factor yheat from equation (3-4). The result is that the direct X-ray heating, in the format of 
equation (2-12), can be expressed as 

= yheat^ionC; (3-5) 

where Tion = ^ion/k. As discussed in Appendix B and §2.4.1, there is also an indirect X-ray heating 
term (Eq. [2-41]) that arises when the thermal effects of collisional excitation and de-excitation of 
the H atom are expressed in terms of Lyman-a cooling (Eq. [B12]). 

The dependence of eion and yhcat on has been studied by many authors, starting with Spitzer 
& Scott (1969) and most recently by Dalgarno, Yan, & Liu (1999), who give an extensive set of 
references to previous work. We have found that the parameterization of Shull & Van Steenburg 
(1983) for atomic H and He mixtures is useful in calculating the effects of X-rays on the inner x-wind 
with moderate mass-loss rates and not too large electron fractions. Their results are confirmed by 
Dalgarno et al. (1999) (who also provide the only theory for situations where the hydrogen is partly 
or fully molecular). According to Shull &; Van Steenburg, the energy to make an ion pair can be 
written, 

6ion /(H) + /(He)' ^^^^ 
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where 

yn = 0.3908 (l - x0-^°92) 1-7592 ^ ^ ^ ^^^^ _ ^o.4614) 1-666 ^ ^3_^^ 

and /(H) and /(He) are the ionization potentials of H and He. The heating fraction is 

„0.2663^ 1-3163" 



yhcat = 0.9971 



1-1 



(3-8) 



It should be noted that existing theories of electron energy loss, on which our ionization and heating 
rates are based, do not hold much beyond Xe = 0.1. In particular, equation (3-5) breaks down in 
the limit — 1, as can be seen from the behavior in this limit of equations (3-7) and (3-8). 
Fortunately, the maximum electron fractions encountered in the present calculations rarely exceed 
Xe = 0.1. 



Putting all of this together for a thermal spectrum with temperature Tx, the ionization rate 
at distance r from an X-ray sources is, 

C«Cx(f^)'(^)/p(rx,eo), (3-9) 
where is the primary ionization rate at a distance r = R^. A useful numerical form for (^^ is 

ATTB^kTx -l-l"^^!" ' l^lOSOergs-iAkevJ \ ) ' ^ ^ 

The function Ip(rx,^o) in equation (3-9) describes the attenuation of the X-rays. As discussed by 
GNI, it decreases more rapidly with optical depth than a power law at large optical depth and 
flattens out to a constant for very small optical depth. In the application to protoplanetary disks, 
the integral was evaluated numerically, but this is infeasible for detailed modeling of the x-wind. In 
Appendix C, wc obtain an asymptotic form for /^(tx, ^o) using the method of steepest descents (as 
did Krolik & Kallmann 1983 and GNI) and then develop a simple interpolation method to combine 
the approximations for small and large optical depths. 



4. Ambipolar Diffusion Heating 

Long familiar from thermal considerations of interstellar clouds (Biermann &: Schliiter 1950; 
Scalo 1977; Mouschovias 1978; Lizano &; Shu 1987), the heating associated with ambipolar diffusion 
(e.g., Mestel & Spitzer 1956; Spitzer 1978; Shu II 1992) was first applied to protostellar winds by 
RGS and to disk winds by Safier (1993). Because RGS worked with a prescribed spherically- 
symmetric wind, they calculated the drag force of the ions on the neutrals from the wind equation 
of motion, i.e., as the net force of gravity plus acceleration. For T-Tauri stars with mass-loss rates 
in the range 10~^ — lO~^M0yr~^, RGS (Figure 9) obtained temperatures in the 4,000 — 5,000 K 
range within 10 but less than 100 K beyond 10^ i?* 10^^ cm (where the temperature decreases 
adiabatically as the 4/3 power of the distance). In this paper, we include ambipolar diffusion 
heating from first principles using the dynamical solution of Shang (1998) described in §2.1. 
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We use an improved approximation for the volumetric rate of ambipolar diffusion heating 
because, unhke the situation in interstellar clouds, the ionization level in the wind may not be very 
small: 

Tao = .. . (4-1) 

Here pn and pi are the mass densities of the neutrals and the ions, respectively, fL is the Lorentz 
force, 

fL = :^(V X B) X B, (4-2) 

and 7 is the ion-neutral momentum transfer coefficient. A short derivation of equation (4-1) is given 
in Appendix D, based on the approximation that the difference in the acceleration (rather than 
the velocity) of the neutrals and ions can be ignored. When pi « p^, equation (4-1) reduces to 
the usual one for low-ionization situations, e.g., Eq. (27.19) of Shu (II 1992). It has the important 
property that it does not become singular as p^ vanishes. 

In Appendix E, wc develop an improved formula for 7 that takes into account the latest 
calculations and experiments on the collision of ions with atomic and molecular hydrogen and 
with helium, including exchange scattering in H"^ + H collisions. For the case of no molecular 
hydrogen and xjie = 0.1, equation (E8) yields 



1Q X 1014 r „,,2 

7 = • " {3.23 + 41.0T^ \l + 1.338 x IQ-^ !^f^}x{B) + 0.243 

1 — 0.714a;e 14 



cm^s-^g-\ (4-3) 



where T4 is the temperature in units of 10,000 K and is the slip speed (w = v; — Vn) in km s^^. 
For temperatures approaching 10^ K, the new 7 is an order of magnitude larger than the value used 
by RGS and elsewhere in the literature. It agrees with Draine's prescription (1980) only for cold 
clouds. 

The slip velocity can be obtained from equation (D8) of Appendix D, 

w = -. (4-4) 

IPAPn + Pi) 

By eliminating 7 from the last two equations, we obtain a quartic equation for w, whose solution 
permits the coupling coefficient 7 to be calculated and then the ambipolar diffusion heating to be 
found from equation (4-1). The contribution from He to 7 in equation (4-3) (the last term) is 
always small. 

A critical factor in the formula for ambipolar diffusion heating is the square of the Lorentz force 
(Eq. [4-2] ) . We calculate the Lorentz force by numerical interpolation on the global x-wind solution 
described in §2.1, which is itself an interpolation between an interior and an exterior (asymptotic) 
solution. Numerical experiments with different interpolation grids indicate that we have accurately 
calculated the Lorentz force for our approximate solution, but we do not have a good estimate of 
the error in the solution itself. However, it is reassuring that our calculations for a model close to 
that of RGS are consistent with their results. 



-20- 



Figure 1 shows the results for what we will call our fiducial model, defined in Table 2, except 
that the parameter a for mechanical heating (to be discussed in §5) has been set equal to zero. In 
other words, Figure 1 includes all of the processes listed in Table 1 except for mechanical heating. 
The upper panels show temperature contours and the lower panels ionization contours in the w-z 
plane. The spatial dimensions are AU and the scale of the plots increases from right to left. On 
the smallest scale (at the extreme right), the range in w is 10 AU and the range in z is 100 AU. 
The temperature very close to the axis approaches 9,000 K, but the region hot enough to excite 
the forbidden lines, roughly T > 6, 000 — 7, 000 K, occupies only a thin inner layer of the jet with 
a thickness less than 1 AU. The extreme right panel corresponds to an angular scale of less than 
0.1" for objects at a distance of 150 pc, one not yet generally available to direct observation. The 
scale has begun to be explored by recent measurements with the Hubble Space Telescope (Bacciotti 
et al. 2000) and with an adaptive optics system on the Canada-Hawaii-France Telescope (Dougados, 
Cabrit, Lavalley, and Menard 2000). On the commonly observed scales shown in the remaining 
panels of Figure 1, ranging from 1/3 to 3 arcsec, the wind is warm at best. The ionization fractions 
range from a few to 10%. 

An analysis of the terms contributing to the basic ionization and heat equations, (2-13) and 
(2-12), reveals that a variety of processes contribute, especially close to the star. For example, 
the gas starting out on the last streamlines closest to the axis is ionized by both X-rays and 
by photoionization of the n = 2 level and the negative ion of the H atom (by stellar and hot-spot 
radiation). The total ionization and recombination rates balance approximately at first, but, within 
a few AU, the ionization falls below the recombination rate and the electron fraction decreases slowly 
with distance (see section §6). For streamlines that start out at larger angles with respect to the 
axis, the X-rays are more important than stellar radiation for ionization. Because the density 
decreases more rapidly with distance along these streamlines, the ionization is essentially frozen 
into the streamline and decreases even more slowly with increasing distance. 

Without mechanical heating. X-rays are generally the most important heating mechanism and 
adiabatic cooling the most important cooling mechanism. However, Lyman-a cooling is larger than 
adiabatic within the inner 10 AU, and ambipolar diffusion heating competes with X-rays at large 
distances beyond 500-1,000 AU. Because the thermal time scale is much shorter than the ionization 
time scale, the general dominance of adiabatic cooling over X-ray and ambipolar heating means that 
the temperature decreases rapidly on every streamline. Of course, as already noted, these heating 
mechanisms are relatively weak. It is perhaps not surprising that X-ray heating is not that effective, 
because only a small fraction of the system luminosity is in X-rays. In accord with equation (4-1), 
ambipolar diffusion heating is weak because the X-rays produce a relatively high level of ionization 
and because the coefficient 7 is an order of magnitude larger than used by previous workers. Much 
higher temperatures, along the lines obtained by Safier (1993) for disk winds, could be achieved 
without X-rays by using the conventional small (but incorrect) 7. For our objective of obtaining 
physical conditions that are compatible with the optical observations of jets, a wind thermal model 
based on stellar radiation. X-rays, and ambipolar diffusion heating is clearly inadequate. 
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5. Mechanical Heating 



The previous section makes clear that neither heating by ambipolar diffusion nor heating by 
X-rays and the other radiation fields in the problem suffice to explain the observed emission lines 
of the abundant heavy elements in YSO winds and jets. We therefore consider whether a small 
fraction of the macroscopic flow energy in the x-wind can be tapped as a volumetric heat source. 
The physical basis for this idea resides in the observation that the actual flows are time-dependent 
(e.g., in the form of pulsed jets as discussed by Raga et al. 1990; Raga and Kofman 1992), with the 
time dependence generating shock waves or turbulent dissipation when fast fluid elements catch 
up with slower ones ^. We will not attempt a detailed discussion of the physics underlying the 
transformation of the fluctuating kinetic energy into heat (which might profitably be studied by 
3-d numerical simulations); instead we appeal to general dimensional reasoning to parameterize the 
functional form of the mechanical heating. 

The volumetric change in the kinetic energy of the flow is represented by the following terms 
in the fluid equations: 



where p and v are the local gas density and flow velocity in an inertial frame at rest with respect 
to the central star. Dimensional analysis suggests that we replace the above expression by 



where s is the distance the fluid element has traveled along a streamline to the location of interest 
in the wind, and where we have introduced a phenomenological coefficient a > to characterize 
the magnitude of the mechanical heating. A choice a <C 1 corresponds to the assumption that 
only a small fraction of the kinetic energy contained in the flow is dissipated into heat via shock 
waves and turbulent decay when integrated over the flow volumes of interest at the characteristic 
distances s in the current problem. 

The expression (5-2) could be made exact if we allowed a to be arbitrarily dependent on the 
spacetime coordinates (x, t) . In practice, we shall make the simplifying assumption that a is a 
global constant, chosen to obtain a reasonable lighting up of the entire wind flow. With a ^ 1, 
the scaling with 1/s in equation (5-2) could then be justified on the basis of the propagation of 
weak planar shocks where the velocity jump across the shock varies asymptotically as the inverse 
square root of the distance traveled (see §95 of Landau & Lifshitz 1959), with the energy deposited 
into heat (in the fixed frame) varying as the square of the velocity jump. In practice, we prefer 



^Hydrodynamic turbulence has often been invoked for heating interstellar clouds (e.g., Black 1987). Rodriguez- 
Fernandez et al. (2001) have recently offered it as an explanation of the warm clouds observed near the Galactic 
Center. McKee and Zweibel (1995) and Ostriker, Gammie, and Stone (1999) discuss the similarities and differences 
between the dissipation of hydrodynamic and MHD turbulence. 




(5-1) 



= a p — . 



(5-2) 
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to regard equation (5-2) not as being derived from specific dissipative processes, but as a generic 
model equation whose form satisfies broad physical considerations and whose utility comes from 
its simplicity of application within the context of small fluctuations about some mean time-steady 
flow. 

From another perspective, when we remember that v tends to v„ ~ 200 — 300 km s^^, it is 
clear that a must be much less than unity, for otherwise the wind will get too hot. We can roughly 
approximate adiabatic cooling by 

^v{-U - -V^- (5-3) 

where r/ is a parameter that is much less than one for streamlines close to the axis and of order 
unity at large angles (77 = 1/3 asymptotically, in the latter case) . When this expression is balanced 
against the mechanical heating equation in (5-2), we obtain a rough estimate for the temperature 

^kT ^ {^)mv\ (5-4) 

From this result, we see that, not only is the collimated jet {rj <C 1) much hotter than the un- 
collimated wide-angle wind, but a must be much less than one in order to avoid heating the wide 
angle wind to extremely high temperatures. The full scale calculations discussed in §6 show that 
a ~ IQ-^ yields jet temperatures in the 5, 000 — 10, 000 K range that have been deduced from 
observations of forbidden lines. 



6. Results 

We now present the results of calculations for a fiducial or reference model, defined in Table 2, 
and for several variations on it. 

Table 2. Fiducial Model 





0.8 Mq 
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4.8 i?. 




7.5 d 




3.2 X 10~* Mq yr 
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4 X 10^^ ergs"-*- 


a 
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As remarked earlier, the parameters have been chosen to represent a solar-mass YSO in a 
fairly active phase. The numerical value of Lx pertains to two sources, one in and one above 
the rcconnection ring, each with a soft and a hard component with individual X-ray powers of 
IqSI ergs~^. The parameter a has been chosen to have the order of magnitude 10^^ on the basis of 
an approximate solution for the temperature that includes only adiabatic cooling and mechanical 
heating and assumes that the density of a collimated streamline varies inversely with the distance. 



6.1. Temperatures and Ionization Fractions 

Figure 2 shows the temperature and ionization profiles (contours of constant T and Xe) on 
various spatial scales in the same way as Figure 1 (for no mechanical heating). It is immediately 
clear from a comparison of the upper panels of Figures 1 and 2 that mechanical heating at this 
level leads to a much warmer wind. The electron fraction in the two sets of lower panels are not 
that different because X-rays dominate the ionization in both cases. The quantitative changes in 
in going from Figure 1 to Figure 2 arise from several temperature-dependent ionization processes: 
recombination (decreases with increasing T), photodetachment of H~ (H~ abundance increases 
with T), and photoionization of the n = 2 level of atomic hydrogen (population increases with T). 

The wide range of physical properties manifested in the x-wind imply that many physical pro- 
cesses are important, although only a few will dominate at any particular location. For the fiducial 
model and modest variations on it, certain processes do play a more global role. For ionization, 
they are X-ray ionization and radiative recombination; for heating, X-rays and mechanical heating; 
and for cooling, Lyman-a and adiabatic cooling. We find that the relative contribution of a process 
varies from streamline to streamline and also with distance along an individual streamline. For 
example, photoionization of the n = 2 level of atomic hydrogen is the most important ionization 
process within a few AU of the star for the last 10% of the streamlines close to the jet axis. X-rays 
then take over and arc effective over distances of 10-20 AU for producing the initial ionization of 
these streamlines. At larger distances, the total ionization rate is smaller than the recombination 
rate. Since the recombination time scale is longer than the dynamical time scale, the ionization 
tends to get frozen into the wind and Xg decreases slowly with increasing distance, as shown in 
the lower left panels of Figures 1-3. This behavior has been found in several jets of young stars, 
e.g., Dougados et al. (2000), Lavalley-Fouquet, Cabrit, and Dougados (2000), Bacciotti, Eisloffel, 
and Ray (1999), and Bacciotti and Eisloffel (1999). For most of the streamlines, the X-rays are 
responsible for setting up the initial ionization of the flow. Similarly, the inner wind out to 10 AU 
is mainly heated by X-rays, but mechanical heating dominates most of the rest of the flow. 

In addition to adiabatic and Lyman-a cooling, gas on the inner streamlines close to the source 
is cooled by several atomic hydrogen processes discussed in §2.4, mainly recombination and H~ 
cooling (equations [2-37] and equations [2-26], respectively). Within a short distance from the 
source, adiabatic and Lyman-a cooling take over and eventually adiabatic cooling dominates. The 
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transition to adiabatic cooling occurs more rapidly for the lower streamlines. The cooling by the 
forbidden lines of 01 and SII also contribute significantly on the inner, collimated streamlines, 
eventually dominating adiabatic for the inner 10% of the streamlines and becoming one-third as 
strong as adiabatic for the inner 25% of the streamlines. For most of the rest of the (uncollimated) 
wind, forbidden line cooling is unimportant. 

Although the wind for the fiducial case (Figure 2) has an ionization fraction in the right range, 
it is not hot enough to emit the forbidden lines at the levels observed in the brightest jets. We 
can achieve the desired temperatures by increasing the coefficient a for mechanical heating. For 
example, Figure 3 shows the temperature and electron fraction profiles when a is increased by a 
factor of two to 2 x 10~^. The temperature is increased by almost a factor of two, and the electron 
fraction by about 20%. 



6.2. Synthetic Images 

Figure 4 shows synthetic images of the wind in the forbidden lines of SII A6731 (left) and 
O I A6300 (right) for the case illustrated in Figure 3, as they might be observed edge-on and with 
near-perfect angular resolution. These images are similar to those constructed earlier by SSG, but 
with some unimportant technical differences in the way that we make the image of the innermost 
regions of the jet. The main physical difference with SSG is that here we calculate the temperature 
and electron fraction from an almost first-principles model, rather than assuming constant values 
for these parameters. The abundance of 01 is calculated on the basis of a theory that gives the 
standard result given by Osterbrock (1989) for HII regions, i.e., O^ is maintained in chemical 
equilibrium by very fast forward and backward charge-exchange (assuming no H2): 

x(0) "g'^P x(H) ^g'^P i-.xe- ^ ^ 

The fact that H+ -|- O charge-exchange is fast means that 0+ comes rapidly into equilibrium with 
the slowly-varying electron fraction. This situation does not generally hold for other atoms where 
charge exchange with H+ is much weaker, especially for sulfur. We are in the midst of developing 
a more general theory of the ionization of the major ionic carriers of jet forbidden lines. In the 
interim, we assume for Figure 4 that all of the sulfur is in SII (as in SSG). 

We obtain the appearance of an optical jet in Figure 4 because the inner x-wind has a stratified 
density profile that varies approximately as the inverse square of vj (SSG), the distance to the jet 
axis, and because the temperature and electron fraction have the right values to produce forbidden 
line emission. The base of the model jet in Figure 4 has a rounded conical shape suggestive of HST 
images (e.g., Eisloffel ei 2000) and a horizontal width of the order of 25 AU. The width, which 
depends somewhat on definition, is determined by the critical density of the transitions, along 
with the temperature and electron fraction. The emissivities integrated along the line of sight 
decline very rapidly with horizontal distance from a cusp close to the axis. A careful inspection of 
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the images near the equatorial plane reveals that 01 A6300 is stronger than SII A6731, basically 
because high-temperature and high-density regions contribute more to the line-of-sight integral of 
the emissivity than the lower temperature and lower density regions that are more important at 
high altitudes. Here we are seeing the result of the competition between the lower critical density 
of the SII A6731 transition and the greater intrinsic strength of the 01 A6300 transition (and 
larger O abundance). Images like Figure 4 provide a concrete basis for observational tests of our 
thermal-chemical theory of the x-wind jet. 



6.3. Additional Parameter-Space Studies 

In addition to the models shown in Figures 2-4, we have made some further exploratory 
calculations without attempting a systematic search of the model parameter space, defined largely 
by mass-loss rate M„, X-ray luminosity Lx, and mechanical heating strength a. For example, we 
have calculated models with larger values of a. Increasing a from 0.002 to 0.005 increases the 
temperature by about 2/3 and the electron fraction by about 1/3. Inside the jet {vj < 25 AU), 
T is in the 10,000 - 13,000 K range and Xg is in the 0.03-0.05 range. This model may well have 
astrophysical applications, and its observational aspects will differ from those of the a = 0.002 
model in Figures 3 and 4. For example, higher-excitation levels of 01 and SII may be excited and 
significant abundances of Oil and Nil produced. An important question is whether a warm model, 
with a = 0.005 or larger, could produce a level of ionization sufficient to produce the forbidden 
lines by just collisional ionization of the H atom without X-rays. Setting Lx = in the model 
with a = 0.005, we find that the electron fraction is reduced by about an order of magnitude and 
that the temperature is increased even further, e.g., up to 15,000 K inside the jet and even higher 
outside. (The decrease in Xe makes ambipolar diffusion heating effective close in and reduces some 
of the cooling). Obtaining wind ionization levels greater than a few percent by heating without 
X-rays doesn't work without going to temperatures beyond the range indicated by the forbidden 
line observations. Some type of external mechanism is required to produce the degree of ionization 
inferred from observations, and we have shown that stellar X-rays are able to do the job. 

When we examine the upper part of Figure 2, we see that the temperature of the wide-angle 
wind is quite high on the largest scale shown (the upper left panel), roughly between 3,000-6,000 K; 
for a = 0.002 (Figure 3) the range is 6,000-9,000 K. We are uncertain of the physical significance 
of this result because there arc no measurements of the wind in this region (at least until larger 
distances are reached where the wind collides with ambient material). One possibility is that 
the mechanical heating formula used so far, equation (5-2), docs not hold for the majority of the 
streamlines which are strongly divergent (where Shockwaves do not propagate even approximately 
according to a planar description). Using other prescriptions (e.g., changing the density dependence 
in equation (5-2) from a linear to a nonlinear dependence), we find that the collimated jet can be 
made warmer and the uncoUimated outer wind colder. 
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The total X-ray luminosity used in the fiducial model (Table 2), Lx = 4 x lO^-^ergs"^, is 
several times larger than the values determined by CHANDRA observations of young solar-mass 
YSOs in Orion (Feigelson 2001, private communication). However, the key physical parameter 
is the X-ray ionization parameter, (^/n, which determines the imprinting of the ionization at the 
base of the wind. Thus the essential model parameter is not Lx alone but something closer to 
Lx/Myf, at least for small values of where X-ray absorption effects are small or moderate. We 
have run models for a = 0.002 (to match Figures 3 and 4) where Lx and are simultaneously 
decreased (and increased) by a factor 3. When Lx and Mw are decreased by 3, the temperature 
is reduced slightly and the ionization fraction increased somewhat more due to the reduction in 
X-ray absorption. The net result is a model which is very similar to the ones shown in Figures 3 
and 4. Increasing Lx and iWw may also work, despite the fact that the electron factor is decreased 
by a factor of 2 because of the increase in X-ray absorption, simply because the emissivity of the 
forbidden lines is determined by the electron density rather than electron fraction. It should be 
recalled that Lx can also be greater than 10^^ ergs~^ in flares and that changes in Lx and My, are 
likely to occur on different time scales. 

The warm and ionized conditions found for the outer wind in Figures 1-3 are in stark contrast 
to the results of RGS because of our inclusion of X-ray ionization and mechanical heating in the 
present calculations. We have already discussed our uncertainty in applying the heating model 
of §5 to the uncollimated flow in the absence of compelling observational information about this 
part of the wind. It is quite possible that the wide-angle wind is not as warm as the above figures 
suggest. Furthermore, our thermal-chemical model needs to be extended in this region to include 
the thermal and chemical effects of molecules. Although not mentioned in §2, we have made a 
preliminary study of molecular hydrogen, mainly to insure that the abundance of H2 is negligible 
in the inner part of the wind. The molecular physics and heavy element chemistry of the wide- 
angle part of the wind are of considerable interest in connection with the detection of H2 jets in 
young embedded YSOs (e.g., Zinnecker, McCaughrean, and Rayner 1998; Stanke, McCaughrean, 
and Zinnecker 1998), and we plan to return to this subject in the near future. 



6.4. Line Ratios 

In addition to synthetic images of the forbidden line emission, like those displayed in Figure 4, 
we can also examine particular line ratios that are sensitive to the underlying physical properties 
of the wind as calculated in this paper. This approach has been widely used for HII regions and 
planetary nebulae (Ostcrbrock 1989), and it has been developed into a diagnostic tool for YSO jets 
by Bacciotti and Eisloffel (1999). Figure 5 is a plot of the S IIA6716/S IIA6731 line ratio vs. the 
SIIA6731/OIA 6300 line ratio based on the same model as the synthetic images in Figure 4, where 
the jet is viewed perpendicular to its axis. The former ratio (the ordinate in the figure) is diagnostic 
of rie and the latter (the abscissa) is sensitive to T. The blue dots are the ratios formed from the 
line intensities for each pixel in the synthetic image of Figure 4, plotted one against the other. The 
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dense concentration of blue points at the top of the "blue cliff" arise from distant wind locations 
with low temperature and electron density, whereas the points at the lower left come from close in 
where temperature and electron density are high. The red asterisks are data for HH objects from 
the compilation of Raga, Bohm, and Canto (1996), and the red circles are data for DG Tau taken 
from Lavalley-Fouquet, Cabrit, and Dougados (2000). The comparison is meant to be illustrative 
and should not be taken too literally. The blue points have been obtained by "lighting up" every 
pixel of a steady-state jet, as modeled by a specific x-wind model with mechanical heating according 
to equation (5-2), and viewing the jet from a single direction. The asterisks are observations of 
a diverse set of HH objects viewed at a variety of angles and with a range of spatial resolution. 
Under these circumstances, we should not expect any more than a general kind of agreement. 

The successful demarcation by the theory of the region in Figure 5 where observed YSO-jet 
line-ratios are found is therefore extremely satisfying. It confirms that the temperature and electron 
density (integrated along the line of sight) of our thermally and mechanically self-consistent x-wind 
model are indeed in the right range to explain the observations of real sources. It is significant 
that this agreement is obtained by adjusting the one free parameter at our disposal, a, since the 
others are constrained by independent observations. We emphasize that we use the same value 
of a in the line-ratio plot of Figure 5 as in the image in Figure 4, a = 0.002. Lavalley-Fouquet 
et al. (2000) have attempted to correlate a similar data set with a complex shock model that employs 
a continuous distribution of shock velocities for each of 5 values of pre-shock density, ranging from 
10^ — 10^ cm~^. The resulting line-ratio diagram then consists of a set of 5 curves or branches, one 
for each pre-shock density. Their analysis does not give the broad range of conditions observed for 
HH objects, but it is more successful with the more uniform set of high spatial-resolution data for 
DG Tau. The latter data appear to support a shock interpretation of the line ratios. 

It is worth commenting that some of the extreme line-ratios, represented by the red circles in 
the lower left part of Figure 5, come from the highest spatial resolution measurements available 
for DG Tau obtained with adaptive optics (Lavalley-Fouquet et al. 2000). It is no coincidence, we 
believe, that the high temperature and electron-density conditions required to produce such ratios 
occur in the theoretical model near the base of the observed flow. By increasing the mass-loss rate 
slightly, these data would lie closer to the main body of the theoretical blue points in Figure 5. The 
observations of Lavalley-Fouquet et al. (2000) reinforce the importance of carrying out spectroscopic 
measurements at a spatial resolution sufficient to probe the extreme conditions close to the source 
of the jet. It is also noteworthy that the three red asterisks in the lower right of the diagram that 
fall outside the envelope of blue dots are all associated with strong bow shocks (Raga et al. 1996). 
The large jumps experienced across strong bow shock are evidently not well represented by the 
simple formula, equation (5-2), with a (weak-shock) value a = 2 x 10~^. It should also not be too 
surprising that some data points lie outside of the theoretical range in Figure 5, considering that 
the blue dots have been calculated for a single viewing angle of one specific jet model, whereas the 
data sample a wide range of objects viewed under different conditions. 

The main purpose of the line-ratio plot in Figure 5 is to ensure that the physical properties of 
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the x-wind, as calculated in this paper, provide a sound foundation for quantitative comparisons 
between theory and observations. More detailed comparisons will require tailoring the theoretical 
models to the specific parameters of individual sources and considering additional diagnostics, e.g., 
the forbidden lines of other species and the radio continuum emission and radio recombination lines 
discussed in the next section. This next level of modeling would better constrain the adjustable 
parameters of the problem as well as test the theory under a wide range of flow and radiation 
conditions. Given the developments of this paper, such detailed tests are within our grasp, but 
their implementation is beyond the scope of the present paper. 



7. Discussion 

The results presented in §6 confirm that the thermal-chemical program described in this paper 
can provide the basis for making comparisons between x-wind theory and observations. Using the 
forbidden lines S II (A6731) and O I (A6300) from jets as the main example, we have shown that the 
collimated portion of the x-wind, when ionized by X-rays and heated mechanically, emits these lines 
in a manner strongly suggestive of the actual images made by observers (Figures 4). The model 
also appears capable of reproducing the measured line ratios (Figure 5). More definitive conclusions 
will require detailed quantitative comparisons between the theory and observations of individual 
objects. Thus we are in the midst of a study that considers further aspects of the forbidden lines, 
such as diagnostic line ratios of additional atoms and ions. The analysis of forbidden line images 
obtained at high spatial and spectral resolution have the potential to provide strong tests of the 
predictions that we are now able to make for the x-wind. In order to explore the jet structure 
predicted for scales smaller than 25 AU, observations with an angular resolution significantly better 
than 0.15" are required for sources at a distance of 150 pc. Recent measurements with the Hubble 
Space Telescope (Bacciotti et al. 2000) and with an adaptive optics system on the Canada-Hawaii- 
France Telescope (Dougados et al. 2000) have begun to probe jets on this scale. 

An independent test of the x-wind model can be made by interferometric measurements of 
the radio continuum emission from the hot partially-ionized gas in the inner wind. Because of 
the absence of extinction, radio observations can probe deep into the inner wind close to the 
source of the outflow. Thermal radio emission has already been detected at the centers of more 
than 100 YSO outflows (Eisloffel et al. 2000; Rodriguez 1997), and about one fifth of these have 
been mapped at high spatial resolution, including a fair number that can be modeled to test the 
x-wind (Rodriguez 2001, private communication). According to the standard theory of thermal 
brcmsstrahlung emission from a plasma with a uniform temperature (e.g., Shu I 1991), the optical 
depth is 

Tff «i 6 X 10-^^Ecm^ (A/3cm)2 T^^''^ , (7-1) 

where is the temperature in 10^ K and E is the emission measure. In the approximation of 
Paper V, where the density at a given height above the midplane varies with horizontal distance 
as ci7~^, Tff oc A^ci7~^, and the size of emission contours of a given intensity level should scale as 
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A^/^, as observed in the best studied cases (see the review of Anglada 1996). Shorter wavelength 
observations are then favored for discriminating the wind from the cooler and lower density inner 
disk. We expect that the thermal brcmsstrahlung emission for models with the fiducial parameters 
in Table 2 will be mainly optically thin, except very close to the source. We plan to synthesize 
the emission for models of the x-wind as it would be observed by the Very Large Array. If we 
are successful in reproducing the observations, we should be able to determine the mass-loss rates 
of thermal jets in a more realistic way than is currently done with the simple bi-conical model of 
Reynolds (1986). 

Another way of testing the x-wind model with radio observations is provided by the mm and 
sub-mm recombination lines emitted by the hot plasma near the base of wind. Important kinematic 
information can be obtained by sensitive measurements of line shapes at high spatial resolution. 
High spectral resolution is also required to detect the lines in the presence of strong continuum 
emission by the disk at mm and sub-mm wavelengths. The hydrogen recombination lines that lie in 
the mm and sub-mm bands for a transitions (n + 1 — ^ n) occur for principal quantum numbers n in 
the range 25-45. Welch and Marr (1987) made the first detection of mm recombination lines in the 
ultra-compact HII region W3(0H) with the H42a line at 86 GHz. This discovery was soon followed 
by other detections in regions of massive star formation (e.g., Gordon and Walmsley 1990) and by 
the discovery of masing transitions in MWC 349 by Martin-Pintado, Bachiller, Thum, and Walmsley 
(1989). Ground based and space observations of MWC 349 show that the masing reaches a broad 
peak at n = 19 or 300 /xm (Thum et al. 1998). To the best of our knowledge, radio recombination 
lines have not yet been detected in T TaTiri stars. Rough preliminary estimates suggest that the 
radio recombination transitions produced in the x-wind will be weakly masing for lines that fall in 
the familiar sub-mm windows. We expect that some of these lines will be detected by new sub-mm 
instrumentation now under development such as the Sub Millimeter Array (SMA) and the Atacama 
Large Millimeter Array (ALMA). In order to calculate the emission and to synthesize images as 
observed with these new instruments, we need to supplement the program described in this paper 
with a full multi-level population calculation for hydrogen and with appropriate radiative transfer. 
These developments are now in progress, and we hope to report soon on the diagnostic prospects 
of the radio recombination lines for testing the validity of the x-wind model. 



8. Conclusion 

We have developed a thermal-chemical program that provides the basis for making detailed 
predictions for the x-wind model that can be compared with observations. The program incorpo- 
rates new physical processes, particularly for heating (mechanical) and ionization (X-rays). The 
rate coefficients for all of the underlying microscopic processes have been re-evaluated and recalcu- 
lated as required. In some cases, significantly different values have been obtained from those in the 
literature, e.g., the coefficient for ambipolar diffusion heating, and these should be useful in other 
problems. 
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In principle the program can describe a wide variety of flows, all within the context of the 
x-wind model. The key astrophysical model parameters are the mass-loss rate (M„), the X-ray 
luminosity (Lx), and the mechanical heating strength (a). The first two can be chosen to represent 
a particular kind of YSO at some stage of evolution, but a considerable range in these parameters 
is allowed by the observations. They may also be variable on short time scales, as is the case for 
the X-ray emission. In contrast, we regard a as a phenomenological parameter. For the case of an 
active but revealed source with an optical jet, the temperature of such jets, as seen in the forbidden 
lines of oxygen and sulfur, indicates that a 2 x 10~^. 

It is very likely that the three parameters, My^, Lx, and a are not all independent of one 
another. For example, we might expect that all three parameters decrease as we proceed from 
very young and active YSOs to older and less active ones. In this paper, we have concentrated on 
sources with optical jets to illustrate our approach. The exploration of other cases should lead to 
new opportunities for testing the x-wind model. In this context, an interesting question is what 
kinds of jets occur at earlier evolutionary stages when the mass-loss rate is much larger than we 
have used here, ~ 3 x lO~^M0yr~^. To answer this question, we are planning to extend the 
underlying physics of our model to include the essential molecular processes that are expected to 
occur when both the X-rays and the stellar radiation are more heavily extincted than in the cases 
treated in this paper. 

The main result of this paper is the demonstration that the x-wind model, when extended to 
include thermal and chemical processes, has the capability to reproduce line ratios as well as images 
of the forbidden lines of jets in a self-consistent manner. We have also outlined a program for testing 
the model against the observations of individual objects, with the goal of better constraining the 
parameters of the problem and testing the theory under as wide a range of conditions as possible. 
The developments undertaken in the present paper now make such detailed tests possible. 

The authors would be pleased to make available digital versions of the theoretically calculated 
emissivities to observers interested in making detailed comparisons between theory and observation. 
This research has been supported in part by the National Science Foundation through collaborative 
research grants to the Berkeley Astronomy and the NYU Physics departments. S.L. acknowledges 
support from DGAPA/UNAM and CONACyT. The authors are grateful to Alex Dalgarno, Dave 
Schultz, and Pcdrag Krstic for advice and help on the cross sections for the interaction of H"*" 
with H, He, and H2, and they would like to to thank Luis Rodriguez for a careful reading of the 
manuscript and for his comments. 
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A. Rate Coefficient for H -H"*" Neutralization 

The exothermic channels of the reaction 

H-+H+^H(1)+H(n) (Al) 

have the energy yields 12.582, 2.648, 0.758, and 0.096 eV for n = 1,2,3 and 4. But curve-crossing 
considerations and detailed theoretical calculations (Fussen &: Kubach 1986) indicate that, at mean 
center of mass energies less than 2 eV, reactions to the n = 3 level dominate by a large margin over 
those to n = 2. 

The total neutralization cross section has been measured from 0.15-300 eV by Moseley, Aberth, 
& Peterson (1970), from 5-2000 eV by Szucs et al. (1984), and from 30-2000 eV by Peart, Bennett, 
Sz Bolder (1985). Only the lowest energies below a few eV are relevant for our astrophysical 
applications, and here the cross sections of Moseley et al. decrease with energy E roughly as E~^. 
However, the later experiments at higher energies clearly show that the results of Moseley et al. are 
too large by a factor of three. This conclusion is supported by the theoretical calculations of Fussen 
&: Kubach (1986) with good potential energy curves. They give an analytic fit to the data below 3 eV 
(after renormalizing the low-energy results), from which we obtain the following approximation to 
the rate coefficient for H~-H+ neutralization valid below 10,000 K: 



2.40 X 10-^r-V2 + 4 96 x 10"*^ + 6.46 x lO^^^T^/^ ^ 7 45 x lO-^^T 



cm'^s"-^ (A2) 



A typical value at T = 10, 000 K is 3.62 x 10 ^cm^s ^, which can be compared to the RGS value, 

IQ-^cm^s^^ 



B. Level Population of the Model Hydrogen Atom 

Wc treat here the effects of coUisional and X-ray ionization and excitation on the population 
of our model two- level (n = 1,2) plus continuum (c) H atom. Table B lists the relevant processes 
and associated rates. 



Table B. H Atom Processes 


CoUisional excitation and de-excitation 


C12, C21 


Spontaneous decay 


^21/321 


Photoionization by the Balmer continuum 


52c 


Radiative recombination 


ks{> 2) 


CoUisional ionization from n = 1, 2 


Clc, C2C 


X-Ray ionization from n = 1, 2 


Clc, C2C 


X-Ray excitation (n = 1 to ra > 1) 


C12 



The coUisional rate coefficients have the form C = Xenk{T) because we assume that electrons 
are the most important collision partners. We have used Voronov's (1997) fit for kic, at T = 10^ K, 
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it agrees to within 10% percent with the simpler formula used by RGS. For k2c, we fit the rate 
coefficients given by Janev et al. (1987) to obtain 

k2c = 7.37 X 10-10rV2e-39,471K/T ^^3 g-1 (gj) 

For T = 10^ K, kic is about 5 x 10~^ smaller than so that, if the n = 2 level is thermally 
populated, collisional ionizations from the n = 2 level proceed at a faster rate than from the 
ground level. We use the same collisional de-excitation rate coefEcient for the n = 2 — 1 transition 
as RGS, based on Vernazza, Avrett, and Loeser (1981). The photoionization rate is given above in 
§2.4.1: g2c = 5A,*(2) + Wh5A,h(2) (Eq. [2-31]), as is the recombination rate, A;s(> 2) (Eq. [2-34]). 

The X-ray rates in Table B are related to the rate ( at which ion pairs are produced, discussed 
in §3. Because the cosmic X-ray absorption cross section in equation (3-1) is normalized to the 
abundance of hydrogen nuclei, the ionization rate per unit volume is ^n, where n is the density 
of total hydrogen (nn), whereas the rates in Table B are defined in terms of level densities. In 
a H/He mixture, 88% of the X-ray produced ions are H+ ions (Dalgarno, Yan, & Liu 1999). If 
we ignore the small contribution from direct X-ray ionization of atoms in the n = 2 level, then 
0.88^n = Cic^^i- Using the approximate conservation relations, n = ni + ric and Ue = ric, leads to 
0.88Cn = Cic?^(l ~ ^e), or 

Cic = (B2) 

Furthermore, for every H+ ion produced, 1.73 excited H atoms are produced (0.37 with n > 2 
(Dalgarno et al. 1999). Taking into account the radiative branching of these excited levels to n = 2, 
the rate Ci2 is given by 

Cl2 = T^C- (B3) 
1- Xe 

The rate equations based on the processes in Table B are: 

(Ci2 + Ci2)ni + XenA;s(> 2)nc = (^21 + A2i/?2i + 52c + C2C + (20)^2 (B4) 

Xenks{> 2)nc = {Cic + Cic)?^i + {92c + C'2c + C2c)n2, (B5) 
ni + n2 + lie = n, ric = Xeii. (B6) 

The ratio of the populations in the n = 1, 2 levels can be obtained by subtracting equation (B5) 
from equation (B4): 

(C12 + Cic + Cic + Ci2)ni = (C21 + A2ip2i)n2, (B7) 

!^^Q= fe2in.Xe + ^21/321 .gg. 

n2 fci2nXe + kicUXe + Cic + Cl2 

The population of the n = 2 level given in equation of §2.1.4 is obtained by ignoring n2 in the 
conservation relation, equation (B6): ni ^ (1 — Xe)n. A quadratic equation for the steady electron 
fraction can also be obtained from equation (B5), but in this work we find Xe by integrating the 
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rate equation (2-13). It may also be noted that recombination into and ionization processes out of 
the n = 2 level drop out of the population ratio, nxjn^ = Q- 

The thermal implications of the above rate equations can be written as, 

r-A = in2C2i-niCi2)E2i+g2cn2^E2c-Xenksi> 2)nckT -CicUiI -C2cn2^E2c+Cic^fini, (B9) 

where Aei is the X-ray heating energy when an ion pair is generated from the n = 1 level, essentially 
the same as the heating energy used in §3, Aehcat = 2/hcat^eion- When we use equation (B2) for 
(^ic, the last term of equation (B9) becomes the usual X-ray heating term (and now includes He+), 
i.e., 

r-A = {n2C2i-niCi2)E2i+g2cn2^E2c-Xenks{> 2)nckT-CicniI-C2cn2AE2c+CnA€hea.t- (BIO) 

We can rearrange this result using population balance equation (B7), 

(n2C2i - niCi2)^2i = -A2i/32in2^2i + (Cic + Cic + Ci2)^2irii, (Bll) 
so as to exhibit the conventional Lyman-a cooling: 

r-A = -^21/521^2^21 + 92cn2AE2c - Xenks{> 2)nckT - ^{CicTii + C2cn2)I 

+ Cn(Aeheat + 2.22^21). (B12) 

We note that this form of the net heating has two X-ray terms, which we may refer to as direct and 
indirect heating. The direct term is the actual X-ray heating, e.g., as calculated by Shull & Van 
Stccnburg (1983) and by Dalgarno, Yan, & Liu (1999), that arises from elastic collisions of (X-ray 
generated) electrons with the atoms and electrons of the partially ionized plasma. The indirect term, 
proportional to £^21 , arises when the net thermal effect of collisional excitation and de-excitation is 
expressed in a form where the conventional Ly-a cooling appears. The large probability for inelastic 
scattering by secondary electron implies that the wind has a diffuse radiation field which we ignore 
for simplicity, except for trapped Ly-a radiation. 



C. X-ray Attenuation 

The integral that appears in equation (3-9), Ip{rx-,^Q)-, describes the attenuation of X-rays for 
a thermal spectrum (Krolik k, Kallmann 1983), 

/•oo 

Ip{rx, ^0) = / exp [- (^ + TxC)] d^, (CI) 

where ^ = E/kTx, Eq is a low-energy cutoff, and the X-ray optical depth is given by 

Tx = (7p^{kTx)N, N = [ ndr. (C2) 

Jo 
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Following Krolik &; Kallman (1983) and GNI, we use the method of steepest descents to derive 
the asymptotic formula 

Ip « f^^^l ^m^"'^/' eM-iU + rx^;nn] ^ Mrx), (C3) 
where satisfies the algebraic equation 

C(P + U=Prx, (C4) 

which we solve using Newton's method with the initial guess Cm = (prx)^/^+^). Derived for large 
(p^Tx + Cm'^)j equation (C3) is accurate as long as > Co, i-e., even for rx as small as 

TO = ^eg(p + Co) = 3.4x10-3 (C5) 

when Co = O-l- For rx < tq, we derive the expansion (valid for Co ^ 1) by repeated integration by 
parts: 

Co H^pW!!^^ l exp[-(Co + rxCo"")] _^., 
^Prom equations (C3) and (C6), we construct the approximate fitting formula for ah rx, 

i'p(Tx,Co) = (^-) i^p(Tx,Co) + (^-) Jpirx). (C7) 



D. Ambipolar Diffusion Heating for Large Ion Fractions 

We analyze ambipolar diffusion heating with a two-fluid model for the ions and neutrals which 
satisfies the equations of motion 

Pnan = Png + fd, (Dl) 
PiBi = Pig - fd + fL, (D2) 

where 

a = (v V)v (D3) 
stands for the steady-flow acceleration for in an inertial frame, g is the gravitational acceleration, 

fd = 7PnPi(vi - Vn) (D4) 

is the volumetric drag force of the ions on the neutrals, and fL = (V x B)/ x B/47r is the volumetric 
Lorcntz force. Addition of equations (Dl) and (D2) yields the equation of motion for the combined 
ion-neutral fluid, 

pa = pg + fL, (D5) 
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where p = pn + Piis the total mass density (without the electrons) and a = {pnSin + PiS^i)/ P is the 
total acceleration. 

The dynamics of the x-wind can be obtained to sufficient accuracy by ignoring the difference 
between and Vi when we use equations (D3) and (D5) to compute the acceleration of the neutrals 
and the ions kinematically (from the mean velocity v = (/JnVn + piVi)/p) and dynamically (from the 
total force pg + fL). We cannot ignore vi — Vn in the heating problem because ambipolar diffusion 
heating vanishes unless the ion and neutral velocities differ. However, we can use the equation 
for the difference in the accelerations obtained by subtracting 1/ pi times equation (D2) from 1/ pn 
times equation (Dl): 

^ + - Vd - -fL- (D6) 



,Pn Pi J Pi 

If we now ignore the difference in the accelerations, but not in the velocities, the left-hand side of 
D6 is zero and this equation yields 

fd = {-^) fL- (D7) 

For lightly ionized media, p\ <^ p^, we obtain fd = fL- This is the expected and familiar result. 
When the ionization fraction is low, the ions have relatively little inertia and quickly reach termi- 
nal velocity governed by the balance of frictional and Lorentz forces, i.e., equation (D2) has the 
approximate solution, = fL. 

We now obtain the slip velocity from equation (D4) for the drag force, 

Vi-Vn= ,\ (D8) 

7Pi(Pn + Pi) 

This generalizes the lightly- ionized expression by replacing pn in the denominator by p^ + p\. 
Expression D8 is asymmetric in n and i because the Lorentz force acts only on charged particles. 
Thus, the slip velocity can become large if p\ becomes very small, but not if p^ becomes very small 
(when the sea of ions, moving under both gravity and the Lorentz force, simply drags the few 
neutrals that are present along with the rest of the almost completely ionized plasma). 

The volumetric rate of heat input into the combined fluid by ambipolar diffusion is 

Tad = fd • (vi - Vn), (D9) 
and substitution of equations (D4) and (D8) leads to 

Tad = (DIO) 
IPiKPu + PiY 

Notice that Fad = (instead of oo) for a completely ionized plasma when pn = but p\ ^ 0. 
Thus equation (DIO) is self-limiting when T approaches and exceeds 10^ K (and the gas becomes 
collisionally ionized. Although the treatment by Safier (1993) differs from ours, his Fad also 
vanishes when pn = 0. 
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E. The Ion-Neutral Coupling Coefficient For Warm HI Regions 

In order to calculate the ambipolar diffusion heating rate from DIO, we need the momentum 
transfer rate coefHcient 7 (dimensions cm^ s~^ g~^) introduced in equation (D4). The velocity field 
for each species consists of a mean velocity v plus a random velocity u, which we assume is thermal. 
Each ion-neutral pair then gives a contribution to 7 which involves a double thermal average of the 
momentum transfer cross section and can be transformed into (Draine 1986): 

fd = IPnPi w = T^jkUjUkm^k ^jkKjk (El) 

where j and k label ionic and neutral species, respectively, w = Vi — Vn, wjk = vj — Vk, rrijk is the 
reduced mass, 

1 

Kjk = —=CjkSj^e^p{-s'jf.) / x^e^p{-x'^)[2xSjkCOsh{2xSjk) -smh{2xsjk)]ajk{xcjk)dx, (E2) 



is a momentum transfer rate coefficient, Cjk = y^2kT/rnjk, and Sjk = w/cjk- We will approximate 
all of the Wjk by the slip speed w of the two-fluid model used in Appendix D and assume that 
the random velocities are all Maxwellian at the same kinetic temperature. When the momentum- 
transfer cross section varies inversely with the velocity, the integral reduces to the usual Langevin 
rate coefficient va{v), whereas for a constant cross section it is (Draine 1986) 



4_ 
-va 



9 ^W^n 

1 H (- 



1/2 

(E3) 



where v = (8kT/nm)^/'^ is the familiar mean relative velocity of a pair with reduced mass m. 
Before evaluating the sum in equation (El), we discuss the cross sections that actually occur in 
this problem. 



E.l. The Cross Sections 

The model described in the main body of the text deals with warm atomic regions where the 
dominant ion is H"*", in contrast to molecular clouds where heavy ions such as C"*" and a variety 
of molecular ions (notably HCO"*") are more important. The main collision partners of the ions 
are H and He atoms. Osterbrock (1961) focused on the central role of the (induced) polarization 
potential, 

in ion-neutral scattering, where ccpoi is the polarizability of the neutral; he calculated the rate coef- 
ficient for momentum transfer between a heavy interstellar-cloud ion and the dominant gas species, 

H, He, and H2. Equation (E4) leads to a temperature-independent rate coefHcient (Langevin- type, 

1/2 

oc a'l)- Following Draine (1980), essentially all authors have taken into account the breakdown 
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of Osterbrock's treatment at high temperatures, due to interactions of shorter range than by 
approximating their contribution by a constant cross section. In the case of Na"*" + H2 scattering, 
for example, Mouschovias & Paleologou (1981) (see also Drainc, Roberge, & Dalgarno 1983) esti- 
mated this as a geometric cross section, 1.67 x 10"^^ cm^. This guess has not been borne out by 
the recent quantum calculations of Flower (2000), where the para-H2-HC0+ cross section varies as 
E~^^^ for energies between 10 — 10^ K, rather than the Langevin E~^/^ dependence, although the 
numerical value of the rate coefficient at 20 K is essentially the Langevin value) . 

Draine (1980) also realized that H+-I-H scattering, basic to our calculations, is strongly affected 
by charge exchange because of the identity of the two nuclei. There are now good experiments of this 
reaction which provide a sound basis for the calculation of the momentum transfer rate coefficient. 
The definitive proof that charge transfer is significant comes from the merged beam experiment by 
Newman et al. (1982), which measures both H"*" + H and H"*" + D scattering with an ion optics 
system that can distinguish between charge exchange and elastic scattering down to energies as low 
as 0.1 eV. They obtain excellent agreement with the theory of Hunter & Kuriyan (1977), which is 
fully quantum mechanical and goes down to 10~^ eV and also agrees with the high-energy theory 
of Dalgarno k. Yadav (1953) and other high-energy experiments (e.g., Gilbody 1994). 

In addition to the definitive work of Hunter &; Kuriyan (1977) for E > 10~^ eV, the momentum 

transfer cross sections have been calculated by Hodges & Breig (1991) in the same energy region, 
and by Krstic & Schultz (1998) for H+ + H, H+ + He, H+ + H2 for energies E > 0.1 eV. It 
is significant that, for none of the basic molecular ions (H2', HeH"*", Hg"), does the scattering 
cross section manifest the pure E~^/'^ dependence of the Langevin theory. But when (diffraction) 
oscillations are averaged out, they all approximate this energy dependence below a certain energy Ei. 
The three cross sections behave differently above Ei, as seen in Table El, which gives approximate 
power law fits to the momentum transfer cross section a^itiE) = auit{Ei){E / Ei)P . The slopes 
below and above Ei are pi and p2, respectively. 



Table E-1. Momentum Transfer Power Law Fits 



System 


El (eV) 


(TmtiEi)^ 


Pi 


P2 


H+ H 


0.01 


165 


-1/2 


-1/8 


H+ + lie 


1.0 


9.8 


-1/2 


-1 


H+ + H2 


5.0 


5.6 


-1/2 


-2 



1. Cross section units: 10 cm^ 



The break points Ei are high enough for H2 and He that the E''^/'^ fit is sufficient for T < 10^ K. 
For H"''-|-H, however, both parts of the fit should be retained, although the high energy or exchange 
scattering dominates for T > 100 K. 
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E.2. Calculations of the Rate Coefficient 



When the momentum-transfer cross sections are approximated by high and low-energy power 
laws, as in Table El, the rate coefficient in equation (E2) consists of two terms, 

K = ai + a2. (E5) 

When we integrate over all energies for the first term, it becomes a constant, Langevin-type, rate 
coefficient, which we tabulate in Table E2. 



Table E2. Rate Coefficients^ 



System 


oti 


H+ + H 


3.23 


H+ + He 


1.52 


H+ + H2 


2.83 



1. Units: 10"^ cm^ s"^ 



For H+ -I- H, we calculate a2(H) by replacing the E~^^^ dependence in Table El by the 
(constant) aAi-crage value for the interval E = 0.01 — 1.0 eV (1.5 x 10~"^'^cm^), integrating over all 
energy, and using Draine's (1986) equation (E3). The result is 

2 

a2 = 4.10 X IQ-^cm^ s"^ Tf^ [1 + 1.338 x 10-^^]^/^ (E6) 

where T4 is the temperature in units of 10,000 K and is the drift speed in km s~^. Notice that, 
for T = 10^ K, a2{H) is more than an order of magnitude larger than ai{H) due to the dominance 
of charge exchange scattering for E > 0.01 eV. 

For the temperature range of interest, 10^ < T, 10^ K, H+ -|- H dominates over H"*" -|- He and 
H+ -I- H2 scattering because of abundance considerations and also because H+ -|- H scattering is 
so much stronger for E > 0.01 eV. The high-energy contribution a2 can be ignored for H2 and He 
because the breakpoints Ei are larger than for H and because the cross sections decrease more 
rapidly than the characteristic low-energy dependence on 1/v. Thus ai in Table E2 gives a good 
approximation (actually upper limit) to the momentum-transfer rate coefficients for H"*" + He and 
H"^ -|- H2 in the temperature region of interest 

On substituting these results into equation (El), we find that the ion-neutral coupling coeffi- 
cient is 

^ ^ 1 x(H)[«(H)2 + a(H)2] + (4/3)x(H2)a(H2)i + (8/5)xHeC«(He)i 
2mH xu + 2x-R2 + '^ne 

or, numerically, 

/2.13 X lO^^N o _i _i 
^ = i 1 - 0.714X. 



2 

{3.23 + 41.0r4°-^(l + 1.338 x lQ-^^f-^}x{B.) + 2.21a;(H2) + 2.43xHe 



(E8) 



-39- 



REFERENCES 

Anglada, G. 1996, in Radio Emission from the Stars and the Sun, eds. A. R. Taylor k, J. M. Paredes, 
(PASP: San Francisco), 3 

Bacciotti, F., Chiudere, C. k Oliva, T. 1995, A&A, 296, 185 

Bacciotti, F. and Eisloffel, J. 1999, A&A, 342, 717 

Bacciotti, F., Mundt, R., Ray, T.P., Eisloffel, J., Solf, J., & Camenzind, M. 2000, ApJ, 357, L49 

Bacciotti, F. 2001, Rev. Mex. AA (Serie dc Confcrcncias), in press 
Biermann, L. k Scliiter, A. 1950, Zs. f. Naturf. 5a, 237 

Black, J.H. 1987, Interstellar Processes, eds. D. J. Hollenbach and H. A. Thronson, Jr., (Dordrecht: 
Reidel), 731-744. 

Calvet, N. k Gullbring, E. 1998, ApJ, 509, 802 

Dalgarno, A., Yan, M., k Liu, W.-H. 1999, ApJS, 125, 237 

Dalgarno, A. k Yadav, H.N. 1953, Proc. Roy. Soc, A66, 173 

Dougados, C, Cabrit, S., Lavalley, C, k Menard, F. 2000, A&A, 357, L61 

Draine, B. T., 1980, ApJ, 241, 1021 

Draine, B. T., 1986, MNRAS, 220, 133 

Draine, B. T., Roberge, W.G., k Dalgarno, A. 1983, ApJ, 264, 485 

Eisloffel, J., Mundt, R., Ray, T. P., Rodriguez, L. F. 2000, in Protostars and Planets IV, eds. V. 
Mannings, A. P. Boss, and S. Russell, (Tucson: University of Arizona Press), 815 

Feigelson E.D. k Montmerle, T. 1999, ARA&A, 37, 363 

Flower, D.R. 2000, MNRAS, 313, L19 

Fussen, D. k Kubach, C. 1986, J. Phys. B, 19, L31 

Garmire, G., Feigelson, E.D., Broos, P., Hillenbrand, L.A., Pravdo, S.H., Townsley, L., k Tsuboi, 
Y. 200, AJ, 120, 1426 

Glassgold, A.E., Mamon, G.A., k Huggins, P.J. 1991, ApJ, 373, 254 

Glassgold, A.E., Najita, J. k Igea, J. 1997 (henceforth GNI), ApJ, 480, 344 

Gilbody, H.B. 1994, Advances in Atomic, Molecular, and Optical Physics, 33,149 



-40- 



Gounelle, M., Shu, F.H., Shang, H., Glassgold, A.E., Rehm, K.E., & Lee, T. 2001, ApJ, 548, 1051 

Gullbring, E., Calvet, N., Muzerolle, J., & Hartmann, L. 2000, ApJ, 544, 927 

Hartigan, P., Bally. J., Rcipurth, B., & Morse, J. A. 2000, in Protostars and Planets IV, eds. V. 
Mannings, A. P. Boss, and S. Russell, (Tucson: University of Arizona Press), 841 

Henke, B.L., Gullikson, E.M., & Davis, J.C. 1993, At. Data Nuc. Data, Tables, 54, 181 

Hodges, X. & Breig, Y. 1991, JGR, 96, 7697 

Hunter, X. k Kuriyan, Y. 1977, Proc. Roy. Soc, A353, 575 
Igea, J. & Glassgold, A.E. 1999, ApJ, 518, 848 

Janev, R. K., Langer, W. D., & Evans, K. 1987, Elementary Processes in H-He Plasmas, (New 
York: Springer) 

Johns-KruU, C. &; Valenti, J.A. 2000, in Stellar Clusters and Associations: Convection, Rotation, 
and Dynamos, eds. R. Pallavicini, G. Micela, and S. Sciortino, (San Francisco: PASP), 371 

Konigl, A. & Ruden, S.P. 1993, in Protostars and Planets III, eds. E.H. Levy and J.I. Lunine 1993, 
(Tucson: University of Arizona Press), 641 

Konigl, A. &; Pudritz, R.E, 2000, in Protostars and Planets IV, eds. V. Mannings, A. P. Boss, and 
S. Russell, (Tucson: University of Arizona Press), 759 

Koyama, K. in Star Formation 1999, ed. T. Nakamoto (Nobeyama: Nobeyama Radio Obs.), 269 

Krstic, P.S. and Schultz, D.R.1998, Atomic and Plasma-Material Data for Fusion 8, 1 

Krolik, J. H. & Kallmann, T. R. 1983, ApJ, 267, 610 

Landau, L. D. & Lifshitz, E. M. 1959, Fluid Mechanics, (London: Pergamon Press) 

Lvalley-Fouquet, C, Cabrit, S., and Dougados, C. 2000, A&A, 336, L41 

Lee, T., Shu, F. H., Shang, H., Glassgold, A. E., Rehm, K. E. 1998, ApJ, 506, 898 

Lizano, S. k Shu, F.H. 1987, in NATO/ASI Physical Processes in Interstellar Clouds, ed. M. Scholer 
(Dordrecht: Reidel), 173 

Gordon, M.A. k Walmsley, CM. 1990, ApJ, 365, 606 

Martm-Pintado, J., Bachiller, R., Thum, C, k Walmsley, C. M. 1989, A&A, 215, L13 
McKee, C.F. k Zweibcl, E. 1995, ApJ, 440, 686 
Mestel, L. k Spitzer, L. Jr., 1956, MNRAS, 116, 50 



-41 - 



Milalas, D. 1970, Stellar Atmospheres (San Francisco: W.H. Freeman) 

Mohanty, S. & Shu, F. H. 2001, ApJ, in preparation 

Morrison, R. & McCammon, D. 1983, ApJ, 270, 119 

Moseley, J., Aberth, W., & Peterson, J.R. 1970, Phys. Rev. Lett.,24, 435 

Mouschovias, T. Ch. 1978, in Protostars and Planets, ed. T. Gehrels, (Tucson: University of Arizona 
Press), 209 

Mouschovias, T. Ch & Paleologou, E. V. 1981, ApJ, 246, 48 
Najita, J, k Shu, F.H. 1994, ApJ, 429, 808 (Paper III) 

Newman, J.H., Cogan, J.D., Ziegler, D.L., Nitz, D.E., Rundel. R.E., Smith, K.A., & Stebbings, 
R.F. 1982, Phys. Rev. A25, 2976 

Osterbrock 1961, ApJ, 134, 270 

Osterbrock 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei, (Mill Valley: Uni- 
versity Science) 

Ostriker, E. C., Gammie, C. F., & Stone, J. M. 1999, ApJ, 513 

Ostriker, E. & Shu, F.H. 1995, ApJ, 447, 813 (Paper IV) 

Peart, B., Bennett, M.A., & Bolder, K. 1985, J. Phys. B18, L439 

Peres, G., Orlanda, S., Reale, F., Rosner, R., k Hudson, H. 2000, ApJ, 528, 537 

Raga, A. C., Canto, J., Binette, L., and Calvet, N. 1990, ApJ, 364, 601 

Raga, A. C. and Kofman, L. 1992, ApJ, 386, 222 

Raga, A. C, Bohm, K.-H., Canto, J. 1996, Rev. Mex. Astron. Astrofis. , 32, 161 

Rodriguez, L. F. 1997, in Hcrbig-Haro Flows and the Birth of Low Mass Stars, lAU Symposium 
No. 182, eds. B. Reipurth and C. Bertout (Dordrecht: Kluwer), 83 

Rodriguez- Fernandez, N.J., Martin-Pintado, J., Fuente, A., de Vincente, P., Wilson, T.L., 
Hiittmeister, S. 2001, A&A, 365, 174 

Reynolds, S.P. 1986, ApJ, 304, 713 

Ruden, S.P., Glassgold, A.E., & Shu, F.H. 1990, ApJ, 361, 546 (henceforth RGS) 
Safier, P. 1993, ApJ, 408, 115 
Seaton, M. J. 1964, MNRAS, 127, 177 



-42- 



Scalo, J. M. 1977, ApJ, 2113, 705 

Schmitt, J. H. M. M., CoUura, A., Sciortino, S., Vaiaba, G. S., Harnden, F. R., Jr, & Rosner, R. 
1990, ApJ, 365, 704 

Szucs, S., Karemera, M., Terao, M., and Brouillard, F. 1984, J. Phys. B17, 1613, 1984 

Shang, H. 1998, Ph. D. thesis, University of Cahfornia at Berkeley 

Shang, H., Shu, F.H., & Glassgold, A.E. 1998, ApJ, 493, L91 

Shu, F. H., Ruden, S. P., Lada, C. J., & Lizano, S. 1991, ApJ, 370, L31 

Shu, F.H. 1991, The Physics of Astrophysics, Vol. I, Radiation, (Mill Valley: University Science 
Books) 

Shu, F.H. 1992, The Physics of Astrophysics, Vol. II, Gas Dynamics, (Mill Valley: University 
Science Books) 

Shu, F.H., Najita, J., Ostriker, E. C., Wilkin, F., Ruden, S., & Lizano, S. 1994a, ApJ, 429, 781 
(Paper I) 

Shu, F.H., Najita, J., Ruden, S., & Lizano, S. 1994b, ApJ, 429, 797 (Paper II) 
Shu, F.H., Najita, J., Ostriker, E. C, & Shang, H. 1995, ApJ, 455, L55 (Paper V) 
Shu, F.H., Shang, H., Glassgold, A. E., & Lee, T 1997, Science, 277, 1475 

Shu, F.H., Najita, J.R., Shang, H., &; Li, Z.-Y. 2000, in Protostars and Planets IV, eds. V. Mannings, 
A.P. Boss, and S. Russell, (Tucson: University of Arizona Press), 789 

Shull, M. & Van Steenburg, M.E. 1983, ApJ, 298, 268 

Skinner, S. L. & Walter, F. M. 1998, ApJ, 509, 761 

Spitzer, L. Jr. & Scott, E. H. 1969, ApJ, 158, 161 

Spitzer, L., Jr. 1978, Physical Processes in the Interstellar Medium, (New York: Wiley-Interscience) 

Standi, P. C, Lcpp, S., & Dalgarno, A. 1998, ApJ, 509, 1 

Stanke, T., McCaughrean, M. J., k Zinnecker, H. 1998, A&A, 332, 307 

Stelzer, B., Neuhauser, R., & Hambaryan, V. 2000, A&A, 356, 949 

Thum, C., Martin-Pintado, J., Quirrenbach, A., & Matthews, H.E., A&A, 333, L63 

Tsuboi, Y. 1999, Ph.D. thesis, Kyoto Univ. 

Vernazza, J.E., Avrett, E.H., & Loeser, R. 1981, ApJS, 45, 635 



-43- 



Verner, E.M., Verner, D.A., Korista, K.T., Ferguson, J.W., Hamann, F., & Ferland, G.J. 1999 
ApJS, 120, 101 

Verner, E.M., Verner, D.A., Balwin, J.A., Ferland, G.J., amd Martin, P.G. 2000, ApJ, 543, 831 

Voronov, G.S. 1997, ADNDT, 65, 1 

Welch, W.J., Jr. & Marr, J. 1987, ApJ, 317, L21 

Wishart, A. W. 1979, MNRAS, 187, 59 

Woitke, P. & Scdlmayr, E., 1999, A&A, 347, 617 

Zinnecker, H., McCaughrean, M.J., k Rayner, J., 1998, Nature, 394, 862 



This preprint was prepared with the AAS IM^^X macros v5.0. 




Fig. 1. — Temperature (upper) and ionization (lower) contours in the for the w — z plane for the 
fiducial case (defined in Table 2), but with a = (no mechanical heating). The units for the spatial 
scales are AU. Note that the spatial scale proceeds from smaller to larger going from right to left. 




Fig. 2. — Temperature (upper) and ionization (lower) contours in the vo — z plane for the fiducial 
case (defined in Table 2). The units for the spatial scales are AU. 




Fig. 3. — Temperature (upper) and ionization (lower) contours in the vo — z plane for the fiducial 
case (defined in Table 2), but with q = 0.002. The units for the spatial scales are AU. 
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Fig. 4. — Synthetic images of the S II A6731 (left) and O I A6300 (right) brightness for the same 
model as in Figure 3 adapting the methods in SSG. The logj^Q of the integrated intensity is plotted 
in units of erg s~^ cui"^ ster~^. 
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Fig. 5.— SII A6716/SII A6731 line ratio vs. the SII A6731/OI A 6300 based on the synthetic 
images in Figure 4, where the jet is viewed perpendicular to its axis. The model parameters are 
given in Table 2 but with a = 0.002. 



